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Abstract 

We show how to perform a fault-tolerant universal quantum computation in 2D architectures 
using only transversal unitary operators and local syndrome measurements. Our approach is 
based on a doubled version of the 2D color code. It enables a transversal implementation of 
all logical gates in the Clifford-|-T basis using the gauge fixing method proposed recently by 
Paetznick and Reichardt. The gauge fixing requires six-qubit parity measurements for Pauli 
operators supported on faces of the honeycomb lattice with two qubits per site. Doubled color 
codes are promising candidates for the experimental demonstration of logical gates since they 
do not require state distillation. Secondly, we propose a Maximum Likelihood algorithm for the 
error correction and gauge fixing tasks that enables a numerical simulation of logical circuits in 
the Clifford-|-r basis. The algorithm can be used in the online regime such that a new error 
syndrome is revealed at each time step. We estimate the average number of logical gates that 
can be implemented reliably for the smallest doubled color code and a toy noise model that 
includes depolarizing memory errors and syndrome measurement errors. 
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1 Introduction 


Recent years have witnessed several major steps towards experimental demonstration of quantum 
error correction [1, 2, 3] giving us hope that a small-scale fault tolerant quantum memory may become 
a reality soon. Quantum memories based on topological stabilizer codes such as the 2D surface code 
are arguably among the most promising candidates since they can tolerate a high level of noise and 
can be realized on a two-dimensional grid of qubits with local parity checks [4, 5, 6]. Logical qubits 
encoded by such codes would be virtually isolated from the environment by means of an active error 
correction and could preserve delicate superpositions of quantum states for extended periods of time. 

Meanwhile, demonstration of a universal set of logical gates required for a fault-tolerant quantum 
computing remains a distant goal. Although the surface code provides a low-overhead implementation 
of logical Clifford gates such as the CNOT or the Hadamard gate [5, 6] , implementation of logical non- 
Clifford gates poses a serious challenge. Non-Clifford gates such as the single-qubit 45° phase shift 
known as the T-gate are required to express interesting quantum algorithms but their operational 
cost in the surface code architecture exceeds the one of Clifford gates by orders of magnitude. This 
large overhead stems from the state distillation subroutines which may require a thousand or more 
physical qubits to realize just a single logical T-gate [7, 8]. Some form of state distillation is used by 
all currently known fault-tolerant protocols based on 2D stabilizer codes. 

The purpose of this paper is to propose an alternative family of quantum codes and fault tolerant 
protocols for 2D architectures where all logical gates are implemented transversally. Recall that a 
logical gate is called transversal if it can be implemented by applying some single-qubit rotations to 
each physical qubit. Transversal gates are highly desirable since they introduce no overhead and do 
not spread errors. Assuming that all qubits are controlled in parallel, a transversal gate takes the 
same time as a single-qubit rotation, which is arguably the best one can hope for. Unfortunately, 
transversal gates have a very limited computational power. A no-go theorem proved by Eastin and 
Knill [9] asserts that a quantum code can have only a finite number of transversal gates which rules 
out universality. In the special case of 2D stabilizer codes a more restrictive version of this theorem 
have been proved asserting that transversal logical gates must belong to the Clifford group [10, 11, 12]. 

To circumvent these no-go theorems we employ the gauge fixing method proposed recently by 
Paetznick and Reichardt [13]. A fault-tolerant protocol based on the gauge fixing method alternates 
between two error correcting codes that provide a transversal implementation of logical Clifford 
gates and the logical T-gate respectively. Thus a computational universality is achieved by combining 
transversal gates of two different codes. A conversion between the codes can be made fault-tolerantly 
if their stabilizer groups have a sufficiently large intersection. This is achieved by properly choosing a 
pattern of parity checks measured at each time step and applying a gauge hxing operator depending 
on the measured syndromes. The latter is responsible both for error correction and for switching 
between two different encodings of the logical qubit. 

Our goal for the first part of the paper (sections 3-7) is to develop effective error models and 
decoding algorithms suitable for simulation of logical circuits in the Clifford-I-T basis. Although a 
transversal implementation of logical T-gates offers a substantial overhead reduction, it poses several 
challenges for the decoding algorithm. First, T-gates introduce correlations between X-type and 
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Z-type errors that cannot be described by the standard stabilizer formalism. This may prevent the 
decoder from using error syndromes measured before application of a T-gate in the error correction 
steps performed afterwards. Second, implementation of T-gates by the gauge hxing method requires 
an online decoder such that a new gauge hxing operator has to be computed and applied prior to each 
logical T-gate. Thus a practical decoder must have running time 0(1) per logical gate independent 
of the total length of the circuit. The present work makes two contributions that partially address 
these challenges. 

First, we generalize the stabilizer formalism commonly used for a numerical simulation of error 
correction to logical Clifford-|-T circuits. Specihcally, we show how to commute Pauli errors through 
a composition of a transversal T-gate and a certain twirling map such that the effective error model 
at each step of the circuit can be described by random Pauli errors, even though the circuit may 
contain many non-Clifford gates. The twirling map has no effect on the logical state since it includes 
only stabilizer operators. We expect that this technique may hnd applications in other contexts. 

Second, we propose a Maximum Likelihood (ML) decoding algorithm for the error correction 
and gauge hxing tasks. The ML decoder applies the Bayes rule to hnd a recovery operator which 
is most likely to succeed in a given task based on the full history of measured syndromes. This 
is achieved by properly taking into account statistics of memory and measurement errors, as well 
as correlations between X-type and Z-type errors introduced by transversal T-gates. Although the 
number of syndromes that the ML decoder has to process scales linearly with the length of the logical 
circuit, the decoder has a constant running time per logical gate which scales as 0(n2”) for a code 
with n physical qubits. The decoder can be used in the online regime for sufficiently small codes, 
which is crucial for the future experimental demonstration of logical gates. A key ingredient of the 
decoder is the fast Walsh-Hadamard transform. A heuristic approximate version of the algorithm 
called a sparse ML decoder is proposed that could be applicable to medium size codes. 

We apply the ML decoder to a particular gauge hxing protocol proposed by Anderson et al [14]. 
The protocol alternates between the 15-qubit Reed-Muller code and the 7-qubit Steane code that 
provide a transversal implementation of the T-gate and Clifford gates respectively. Numerical simu¬ 
lations are performed for a phenomenological error model that consists of depolarizing memory errors 
and syndrome measurement errors with some rate p. Following ideas of a randomized benchmark¬ 
ing [15, 16] we choose a Clifford-|-T circuit at random, such that each Clifford gate is drawn from 
the uniform distribution on the single-qubit Clifford group. The circuit alternates between Clifford 
and T gates. The quantity we are interested in is a logical error rate dehned as pi = I/s', where 
g is the average number of logical gates implemented before the hrst failure in the error correction 
or gauge hxing subroutines. Here g includes both Clihord and T gates. The sparse ML decoder 
enables a numerical simulation of circuits with more than 10,000 logical gates. For small error rates 
we observed a scaling pi = Cp^ with C ~ 182. Assuming that a physical Clihord-I-T circuit has an 
error probability p per gate, the logical circuit becomes more reliable than the physical one provided 
that pl < Pi that is, p < po = C~^ ~ 0.55%. This value can be viewed as an “error threshold” 
of the proposed protocol. The observed threshold is comparable with the one calculated by Brown 
et al [17] for a gauge hxing protocol based on the 3D color codes which were recently proposed by 
Bombin [18]. We note however that Ref. [17] studied only the storage of a logical qubit (no logical 
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gates). We anticipate that the tools developed in the present paper could be used to simulate logical 
Clifford+T circuits based on the 3D color codes as well. It should be pointed out that the threshold 
Pq = 0.55% is almost one order of magnitude smaller than the one of the 2D surface code for the 
analogous error model [19]. This is the price one has to pay for the low-overhead implementation of 
all logical gates. 

The numerical results obtained for the 15-qubit code call for more general code constructions 
that could achieve a more favorable scaling of the logical error rate. In the second part of the paper 
(sections 8-12) we propose an inhnite family of 2D quantum codes with a diverging code distance that 
enable implementation of Clifford-I-T circuits by the gauge hxing method. The number of physical 
qubits required to achieve a code distance d = 2t + l is n = 2t^ + 8t‘^ + 6t — 1. For comparison, the 3D 
color codes of Ref. [18] require n = + 6t‘^ + 4:t + l physical qubits. The new codes can be embedded 

into the 2D honeycomb lattice with two qubits per site such that all syndrome measurements required 
for error correction and gauge hxing are spatially local. More precisely, any check operator measured 
in the protocol acts on at most six qubits located on some face of the lattice. As was pointed out 
in Ref. [18], the gauge hxing method circumvents no-go theorems proved for transversal non-Clihord 
gates in the 2D geometry [10, 11, 12] since the decoder that controls all quantum operation may 
perform a non-local classical processing. 

The key ingredient of our approach is a doubling transformation from the classical coding theory 
originally proposed by Betsumiya and Munemasa [20] . Its quantum analogue can be used to construct 
high-distance codes with a special symmetry required for transversality of logical T-gates. Namely, 
a quantum code of CSS type [21, 22] is said to be triply even (doubly even) if the weight of any 
X-type stabilizer is a multiple of eight (multiple of four). Any triply even CSS code of odd length 
is known to have a transversal T-gate. The doubling transformation combines a triply even code 
with distance d — 2 and two copies of a doubly even code with distance d to produce a triply 
even code with distance d. Our construction recursively applies the doubling transformation to the 
family of regular color codes on the honeycomb lattice [23] such that each recursion level increases 
the code distance by two. The regular color codes are known to be doubly even [24] (in a certain 
generalized sense). Producing a distance-d triply even code requires {d — l)/2 recursion levels that 
combine color codes with distance 3,5,... ,d. We refer to the new family of codes as doubled color 
code since the construction relies on taking two copies of the regular color codes. It should not 
be confused with the quantum double construction from the topological quantum held theory [25]. 
Doubled color codes have almost all properties required for implementation of logical Clifford-|-T 
circuits by the gauge hxing method. Namely, a doubled color code has a transversal T-gate and can 
be converted fault-tolerantly to the regular color code which is known to have transversal Clihord 
gates [24]. Unfortunately, the doubling transformation does not preserve spatial locality of check 
operators. Even worse, some check operators of a doubled color code have very large weight and 
their syndromes cannot be measured in a fault-tolerant fashion. Converting the doubled color codes 
into a local form is our main technical contribution. This is achieved in two steps. First we show 
how to implement all levels of the recursive doubling transformation on the honeycomb lattice with 
two qubits per site such that almost all check operators of the output code are spatially local. We 
show that each of the remaining non-local checks can be decomposed into a product of local ones 
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by introducing several ancillary qubits and extending the code properly to the ancillary qubits. 
This technique is reminiscent of perturbation theory gadgets that are used to generate effective low- 
energy Hamiltonians with long-range many-body interactions starting from a simpler high-energy 
Hamiltonian with short-range two-body interactions [26]. 

The 15-qubit code studied in Section 7 can be viewed as the smallest example of a doubled color 
code. Furthermore, the 49-qubit triply even code with distance d = 5 discovered by an exhaustive 
numerical search in Ref. [27] can be viewed as a doubled color code obtained from the regular color 
code [[17,1,5]] on the square-octagon lattice via the doubling transformation. The 49-qubit code is 
optimal in the sense that no distance-5 code with less than 49 qubits can be triply even [20, 27]. 
Thus the family of doubled color codes includes the best known examples of triply even codes with 
a small distance. 

Although this paper focuses on codes with a single logical qubit, our fault-tolerant protocols can 
be incorporated into the lattice surgery method based on the regular color codes [28]. The former 
would provide implementation of logical single-qubit rotations decomposed into a product of Clifford 
and T-gates while the latter enables logical CNOTs and can serve as a quantum memory. We note 
that efficient and nearly optimal algorithms for decomposing single-qubit rotations into a product of 
Clifford and T-gates have been proposed recently [29, 30]. 

To make the paper self-contained, we provide all necessary background on quantum codes of 
CSS type, the gauge hxing method, and transversal logical gates in Sections 2-4. The effective error 
model describing the action of transversal T-gates on random Pauli errors is developed in Section 5. 
We describe the ML decoding algorithm suitable for simulation of Clifford-|-T circuits in Section 6. 
Numerical simulation of random Clifford-I-T logical circuits based on the family of 15-qubit codes is 
described in Section 7. This section also serves as an example illustrating the general construction 
of doubled color codes. The latter is described in Sections 8-12. First, we highlight main properties 
of the doubled color codes in Section 8. Dehnition of regular 2D color codes and their properties 
are summarized in Section 9. Doubled color codes and their embedding into the honeycomb lattice 
are dehned in Sections 10,11. The most technical part of the paper is Section 12 explaining how to 
convert doubled color codes into a spatially local form. 

2 Notations 

Let F 2 be the u-dimensional linear space over the binary held F 2 . A vector x G F 2 is regarded as 
a column vector with components Xi,... ,x„. We shall write x'^ for the corresponding row vector. 
The set of integers {1, 2,... ,?7,} will be denoted [u]. Let supp(x) C [n] be the support of x, that is, 
the subset of indexes i such that x* = 1. We shall often identify a vector x and the subset supp(x). 
Let |x| = |supp(x)| be the weight of x, that is, the number of non-zero components. Given a subset 
A C [n] let xa G F^"^' be a restriction of x onto A, that is, a vector obtained from x by deleting all 
components Xi with i ^ A. Conversely, given a vector y G F^"^' let y[A] G F 2 be a vector obtained 
from y by inserting zero components for all i ^ A. For example, if n = 5, A = {1, 3, 5}, and y = (111) 
then y[A\ = (10101). The space F 2 is equipped with a standard basis e^,e^,... ,6", where = l[j] 
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is the vector with a single non-zero component indexed by j. We shall use notations 0 and 1 for the 
all-zeros and the all-ones vectors. The inner product between vectors x,y E ¥2 is dehned as 


x^y = '^Xiyi (mod 2). 

i=l 

A linear subspace spanned by vectors x^,, x™' G F 2 will be denoted {x^,..., x^). Given a subset 
A C [n] and a subspace S C let iS[A] C F 2 be the subspace spanned by vectors y[A] with 
y E S. Let F” and be the subspaces of F 2 spanned by all even-weight and all odd-weight vectors 
respectively. We shall use shorthand notations F = F” and O = O” whenever the value of n is clear 
from the context. Given a linear subspace 5 C F 2 let iS-*- be the orthogonal subspace, 

iS-*- = {x G F 2 : x'^y = 0 for all y G 5} 

and <S be the subspace spanned by all even-weight vectors orthogonal to S, 

s = s^ns. 

A subspace S is self-orthogonal if 5 C iS-*-. We shall use identities (*5-'-)-'- = S, (5 -|- T)-^ = S-‘- nT-‘-, 
and = (1). Given a subspace S, let d(S) be the minimum weight of odd-weight vectors in 

d(S) = min {|/| : f eS^HO}. (1) 

Gonsider now a system of n qubits and let Xj, Yj, Zj be the Pauli operators acting on a qubit j 
tensored with the identity on the remaining qubits. Given a vector / G F 2 and a single-qubit Pauli 
operator P let P(/) be the n-qubit operator that applies P to each qubit in the support of /, 

p(f)= n 

jesupp(f) 

For any subset 5 C F 2 dehne the corresponding set of Pauli operators P(S) = {P(/) : / G iS}. A 
pair of subspaces P C F 2 dehnes a group of n-qubit Pauli operators 

CSS(A,B) = (X(A),Z(B)}. 

Any element of GSS(^, P) has a form i'^X{f)Z{g) for some f E A, g E B, and some integer m. 
Note that the group GSS(^, P) is abelian iff A and P are mutually orthogonal, A C since 

Given a vector / G F 2 , let |/) = |/i ® ® fn) be the corresponding basis state of n qubits. Note 

that |/)=X(/)|b). 
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3 Subsystem quantum codes and gauge fixing 

This section summarizes some known facts concerning subsystem quantum codes of Calderbank- 
Steane-Shor (CSS) type [21, 22] and the gauge fixing method. A quantum CSS code is constructed 
from a pair of linear subspaces C F 2 that are mutually orthogonal, . Such pair defines 

an abelian group of Pauli operators S = CSS(^, i3) called a stabilizer group. Elements of S are 
called stabilizers. We shall often identify a CSS code and its stabilizer group. A subspace spanned 
by n-qubit states invariant under the action of any stabilizer is called a codespace. A projector onto 
the codespace can be written as 

n=jiiEG. (2) 

' ' Ge-s 

In this paper we only consider a restricted class of CSS codes such that all vectors in A and B have 
even weight whereas the number of physical qubits n is odd, 

n = 1 (mod 2), A^ 8 ., B C 8. (3) 

We shall only consider codes with a single logical (encoded) qubit. Operators acting on the encoded 
qubit are expressed in terms of logical Pauli operators 

Xl = X{1\ W = y(T), Zl = Z{1). (4) 

Logical operators commute with any stabilizer due to Eq. (3) and thus preserve the codespace. 
Furthermore, Xi,Yl, Z^ obey the same commutation rules as single-qubit Pauli operators X,Y,Z 
respectively. More generally, a Pauli operator P is called a logical operator iff it coincides with W, 
or Zl modulo stabilizers, that is, P = X(a)Z(b), where a G A+al, b G B+(31, and at least one of the 
coefficients a, (3 is non-zero. A logical state encoding a single-qubit state rj = {1 /2){I +aX+(dYP'yZ) 
is defined as 

PL = Pl{v) = OlH, Ol = S{I + + /3Yl + jZl), (5) 

where 5 is a coefficient responsible for normalization Tt{pl) = 1. 

Pauli operators that commute with both stabilizers and logical operators generate a group Q 
called a gauge group. Elements of Q are called gauge operators. Note that a Pauli operator X(/) 
commutes with all stabilizers iff / G B^ . Likewise, X{f) commutes with the logical operators iff 
1^/ = 0, that is, f E 8. Thus X{f) is a gauge operator iff / G i3, where we use notations of 
Section 2. The same reasoning shows that Z{g) is a gauge operator iff G M. Thus the gauge group 
corresponding to a stabilizer group S = CSS(M, 13) is given by 

g = CSS{B,A). (6) 

By definition, S Y Q. Note that EpiE^ = pi for any E E G, that is, gauge operators have no effect 
on the logical state. 

A noise maps the logical state to a probabilistic mixture of states , where E^ are some 

n-qubit Pauli operators called memory errors or simply errors. We shall only consider noise that 


8 


can be described by random Panli errors. Note that errors and Ep have the same action on any 
logical state whenever E^Ea G Q, that is, gange-eqnivalent errors can be identihed. We shall say that 
an error E is non-trivial ii E ^ Q. Errors are diagnosed by measnring eigenvalnes of some stabilizers. 
An eigenvalne measnrement of a stabilizer X{f) has an ontcome (—where ^{f) G F 2 is called 
a syndrome of X{f). A corrnpted state Ep^E"^ with a memory error E = X{a)Z{b) has a syndrome 
^{f) = If reveals whether E commntes or anti-commntes with X{f). A syndrome of a stabilizer 
Z{g) is dehned as ({g) = g'^a. It reveals whether E commntes or anti-commutes with Z{g). By 
dehnition, gauge-equivalent errors have the same syndromes. 

In some cases stabilizer syndromes cannot be measured directly (for example, if a stabilizer has 
too large weight) but they can be inferred by measuring eigenvalues of some gauge operators. For 
example, suppose a stabilizer X(/) can be represented as a product of X-type gauge operators X{g^), 
that is, f = g^ + g"^ + ...+ g"^ for some f & A and g°‘ G B. An eigenvalue measurement of the gauge 
operator X{g°‘) has an outcome (—where ^{g°‘) is called a gauge syndrome. Once all gauge 
syndromes have been measured, the syndrome of X(/) is inferred from ,^(/) = ■^(5'^)+- • •+■^(5'™')- 
Since gauge operators commute with stabilizers, these measurements do not affect syndromes of any 
Z-type stabilizers. However, measuring gauge syndromes may change the encoded state and the 
latter no longer has a form E^pnEa (since some gauge degrees of freedom have been hxed). We shall 
assume that each syndrome measurement is followed by a twirling map 

Geg 

that applies a random gauge operator G drawn from the uniform distribution on Q. The twirling map 
restores the original form of the encoded state E^PlE},^ by bringing all gauge degrees of freedom to 
the maximally mixed state. Gauge syndromes C(fi'“) for ^-fype gauge operators g°‘ ^ A are dehned 
analogously. 

A memory error E is said to be undetectable if it commutes with any element of S. Equivalently, 
E has zero syndrome for any stabilizer. Errors E that are both undetectable and non-trivial [E ^ Q) 
should be avoided since they can alter the logical state. A code distance d is dehned as the minimum 
weight of an undetectable non-trivial error. The code CSS(^, H) has distance 

d = vam {d{A), d{B)}, (7) 

where d{A) and d{B) are dehned by Eq. (1). 

Next let us discuss a gauge fixing operation that extends the stabilizer group to a larger group or, 
equivalently, reduces the gauge group to a smaller subgroup. Consider a pair of codes with stabilizer 
groups 5 = CSS(^, B) and 5' = CSS(^', B') such that ^ C and H C B'. Let Q = CSS(H, A) and 
Q' = CSS(H',^') be the respective gauge groups. We note that S C S' and Q ^ Q'. Let us represent 
^ as a disjoint union of cosets GaG' for some hxed set of coset representatives Gi, . .. ,Gm G Q. 
Simple algebra shows that the codespace projectors H and H' of the two codes are related as 

m 

n = ^ GJVGl (8) 

a=l 
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and all subspaces GJ^G^ are pairwise orthogonal. Substituting Eq. ( 8 ) into Eq. (5) shows that any 
logical state pi of the code S can be written as 

-i m 

Pi = - E P) 

ra 

a=l 

where is a logical state of the code S' . Furthermore, pi and encode the same state. Since 
Ga are not gauge operators for the code S', they should be treated as memory errors. Moreover, 
any pair of errors Ga and Gb can be distinguished by measuring a complete syndrome of the code S' 
(that is, measuring syndromes of some complete set of generators of S'). Indeed, if Ga and Gb would 
have the same syndromes for some a ^ b then the operator GaGl would commute with any element 
of S' as well as with the logical operators. This would imply GaGl G Q' which is a contradiction 
since GaG' 7 ^ GbG'■ Thus a complete syndrome measurement for the code S' projects the state pl 
in Eq. (9) onto some particular term Gap'i^l and the coset GaG' can be inferred from the measured 
syndrome. Applying a recovery operator R chosen as any representative of the coset GaG' yields 
a state RGaP'i,GlB) = p'^. This completes the gauge hxing operation. The gauge hxing can be 
performed in the reverse direction as well. Namely, the encoded state p'^ can be mapped to pi by 
applying randomly one of the operator Gi,..., Gm drawn from the uniform distribution, see Eq. (9). 
We shall collectively refer to the gauge hxing and its reverse as code deformations. 

It will be convenient to distinguish between subsystem and regular CSS codes (a.k.a. subspace or 
stabilizer codes). The latter is a special class of CSS codes such that the stabilizer and the gauge 
groups are the same. In other words, a stabilizer group CSS(M, i3) dehnes a regular CSS code iff 
A = B and B = A. Note that applying the dot operation twice gives the original subspace, A = A, 
whenever AG S and n is odd. Thus A = B iO. B = A. A regular CSS code has a two-dimensional 
codespace with an orthonormal basis 

|0i) = |.4|-‘-'^EI/> |lL>=W|0i> = |.4r‘/^El/ + i)' (“) 

f&A f&A 

such that n = |0L)(0i| -|- |li)(l 2 ,|. In the case of subsystem CSS codes, the gauge group is strictly 
larger than the stabilizer group, so that the codespace can be further decomposed into a tensor 
product of a logical subsystem and a gauge subsystem. 

4 Transversal logical gates 

In this section we state sufficient conditions under which a quantum code of CSS type has transversal 
logical gates. Recall that a single-qubit unitary gate V is called transversal if there exist a product 
unitary operator Uaii = Ui ® U 2 ® ® Un that preserves the codespace and the action of Uaii on 

the logical qubit coincides with V. More precisely, let Pl(p) be a logical state that encodes a single¬ 
qubit state p, see Eq. (5). Then we require that Uaii PLiji)UGii ~ PL(yv'G~^) for all p. To simplify 
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notations, we shall write Uaii = 11^=1 meaning that Uj acts on the j-th qubit. We shall use the 
following set of single-qubit logical gates: 


H 


72 [1 





1 0 
0 i 


and T 


1 0 
0 


It is known that H and S generate the full Clifford group on one qubit, whereas H, S', T is a universal 
set generating a dense subgroup of the unitary group. The set H, S, T is known as the Clifford-I-T 
basis. To state sufficient conditions for transversality we shall need notions of a doubly-even and a 
triply-even subspace [20]. 

Definition 1. A subspace A ^¥2 is doubly-even iff there exist disjoint subsets C [n] such that 


1/n M"*"! — 1/n M 1 = 0 (mod 4) for all f e A. 


( 11 ) 


Definition 2. A subspace A ^¥2 is triply-even iff there exist disjoint subsets C [n] such that 


\f n M~^\ — \f n M 1 = 0 (mods) for all f e A. 


( 12 ) 


We require that at least one of the subsets is non-empty, since otherwise the dehnitions 
are meaningless. Below we implicitly assume that each doubly- or triply-even subspace is equipped 
with a pair of subsets satisfying Eq. (11) or Eq. (12) respectively. In this case A is said to be 
doubly- or triply-even with respect to M^. The original dehnition of triply-even codes given in [20] 
is recovered by choosing = [n] and M~ = 0. 

Consider a pair of disjoint subsets C [n] such that 

m = \M^\ — \M~ \ = 1 (mod 2) 


and dehne n-qubit product operators 

n 

Tau= n n n n h,ii=i[hj. (i3) 

jeM+ j&M- j&M+ j&M- j=l 

The following lemma provides sufficient conditions for transversality of T, S, and H gates. 

Lemma 1. Consider a quantum code CSS(^, B), where S C F 2 are mutually orthogonal subspaces 
satisfying Eq. (3). The code has a transversal gate T™, 

T^uPL{v)Tfi! = PL{T^vT-n, (14) 

whenever B = A and A is triply even with respect to M^. The code has a transversal gate S'™, 

SaaPL{v)Sfil = PL{S^vS-n, (15) 

whenever AffB and A is doubly even with respect to . The code has a transversal H-gate, 


Hall pLijl)Haii = pL{Hr]H), (16) 

whenever A = B. 
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Proof. We start from Eq. (14). Since B = A, the codespace is two-dimensional with an orthonormal 
basis |0 l), |1l) dehned in Eq. (10). Clearly, Tau\f) = e*’^/'^(hnM+|-|/nM |)|j^ state |/) 

of n-qubits. Combining this identity and dehnition of the logical state lO^,) one gets 

TaiilOL) =T,uJ2\f) = = |0i). 

/eyl f&A 

Here we ignored the normalization of |0i). Using the identity |(/-|- 1) nM^| = \M^\ — |/nM^| and 
dehnition of the logical state |1 l) one gets 

Ta«|U) = ^ 1/+ T) = + i) = 

/e^ f&A 

This proves Eq. (14). Consider now Eq. (15). Let PL = CSS(^,H) be the stabilizer group. First we 
claim that PL is invariant under the conjugated action of Sail, that is, 

S..„HSJ=H. ( 17 ) 

Indeed, since Sail commutes with Z-type stabilizers, it suffices to check that SaiiX{f)Spil G PL for all 
f ^ A. Using the identities SXS~^ = iXZ and S~^XS = —iXZ one gets 

S.uX(f)S-,l = il7n«*l-l7™-|A'(/)Z(/) = A-(/)Z(/) e H. 

Here we used the assumption that A is doubly even and A ^ B which implies Z{f) E PL for any 
f E A. This proves Eq. (17). We conclude that Sail preserves the codespace, Saii^Sfij = H. 
Therefore 

Sail XlUS^^I = Sail XlSJU = z^XlZlII and Sail Zi^Sfil = Z^H. 

The assumption that m is odd implies S"^XS~"^ = i^XZ and S'^ZS~'^ = Z. We conclude that Sail 
implements a gate S'™ on the logical qubit which proves Eq. (15). 

The same arguments as above show that PfaiiAHaii = H. Since Haii interchanges the logical 
operators Xl and Z^, this proves Eq. (16). □ 

We note that T = T^S^Z"^ whenever the integers p, q, r obey p -|- 2g -|- 4r = 1 (mod 8). Using this 
identity one can implement the logical T-gate as a composition of the gate with any odd p, the 
gate S'^ with any odd q, and the logical Pauli operator Z^. 

5 Commuting Pauli errors through T-gates 

Consider a regular code CSS(^, H) with the logical states |0 l), |1l) and suppose the code has a 
transversal T-gate that can be implemented by a product operator Tq;; = Let 1-0) = a\Qi)+(3\li) 
be some logical state and p = |'0)('0|- Iii fhe absence of errors Tau preserves the codespace, so that 
r] = TaiipThi is also a logical state. Consider now a memory error X{e). Our goal is understand what 
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happens when Tail is applied to a corrupted state X(e)pX(e). We will show that a composition of Tail 
and a certain Pauli twirling map transforms the initial state X{e)pX{e) to a probabilistic mixture 
of states ErjE\ where E = X{e)Z{f) and / C e is random vector drawn from a suitable probability 
distribution P(/|e). However, this is true only if the initial error e satishes a technical condition 
that we call cleanahility. To state this condition, represent the binary space F2 as a disjoint union of 
cosets of A. By dehnition, each coset is a set of vectors e + ^ for some e G F2. Since p is stabilized 
by X{A), the state X{e)pX{e) depends only on the coset e + A. 

Definition 3. A coset of A is cleanable iff it has a representative e such that no vector g G A-^ fl O 
has support inside e. 

Recall that the set A^ fl O describes undetectable non-trivial Z-errors. Thus a coset is cleanable 
iff it has a representative whose support contains no such errors. In particular, a coset is cleanable 
whenever it has a representative of weight less than d{A), see Eq. ( 1 ). In practice, one can create 
a list of all cleanable cosets by examining all vectors e G F 2 and checking whether a set M(e) = 
{g G : g[e] G A'^} is non-empty. The set M(e) is determined by a linear system over F 2 with 
|e| variables and 1 -|- dim (^) equations. If M(e) is empty, the coset e -1- ^ is marked as cleanable. 
Cosets remaining unmarked at the end of this process are not cleanable. 

From now on we assume that the coset e ^ is cleanable and e is a representative that obeys 
the condition of Dehnition 3. Dehne a twirling map W_a that applies a random X-stabilizer drawn 
from the uniform distribution, 

geA 

Note that IT 4 has trivial action on any logical state. Consider states 

V = TaiipTlii and p = WA{TaiiX{e)pX{e)Tlif). 

The identity TXTl = implies TaiiX{e)Tlii ~ X{e)S'^{e), where ~ stands for some phase 

factor. Thus 

p = WAiX{e)S{eypS{e)X{e)). ( 18 ) 

Next, the identity S (J — iZ) /-\/2 implies 

/Ce 

where / C e is a shorthand for supp(/) C supp(e). This yields 

X{e)pX{e) = 2-1^=' 

fJ'Qe 

We note that p is invariant under all stabilizers X{g) that appear in the twirling map since it is a 
logical state. Commuting a stabilizer X{g) from the twirling map towards p gives an extra phase 


13 


factor Summing up this phase factor over g & A gives a non-zero contribution only if 

f + f e A-^. This shows that 

X(e)77X(e) = 2-1^1 Zif)gZin. (19) 

fJ'Qe 

By cleanability assumption, the sum in Eq. (19) gets non-zero contributions only from the terms 

with f + f' E A'^ E] £ = A. Furthermore, A = B since we assumed that the code is regular. Dehne 
a subspace 

Bie) = {geB : gCe}. 

Performing a change of variables /' = / -|- in Eq. (19) gives 

X{e)fiX(e) = 2-<’<J2 E 

/Ce g&B(e) 

Here we noted that gZ{f + g) = gZ{f) since Z{g) is a stabilizer and ?7 is a logical state. Next, the 
identity \f + g\ = \ f\ + \g\ — 2\f fl g\ implies = (—i)P9+l9l/2_ Note that \g\ is even since 

B C £. We conclude that 

/Ce g,hGB{e) 

Here we introduced a dummy variable h E B{e) and used the fact that rj is invariant under stabilizers 
Z{h) with h G B{e). The sum over h gives a non-zero contribution only if (? € B{e)-^. We arrive at 

A(e)f,A'(e) = ^ P(/|e) Z(/),,Z(/). (20) 

/Ce 

where 

F(/|e) = 2 -l"l (21) 

g£B{e)nB{e)-^ 

We claim that P(/|e) is a normalized probability distribution on the set of vectors / C e, that is, 

5^F(/|e) = l and P(/|e) > 0. 

/Ce 

Indeed, normalization of P(/|e) follows trivially from Eq. (21). To prove that P(/|e) > 0 choose an 
arbitrary basis P(e) fl P(e)“*“ = (^f^,..., g^). For any vector g = one has 

m m 

Ifi'l = - 2 ^ XaXb\g°‘r\g’'\ = ^Xa\g''\ (mod 4) 

a=l l<a<b<m a=l 
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since all overlaps \g°' fl g^\ must be even. It follows that 

m 

P{f\e) = 2-1^1 n + (_l)r9“+l9“l/2] > 0. (22) 

a=l 

To conclude, the transversal operator Tail followed by the twirling map transforms the initial 
error X{e) to a probabilistic mixture of errors X{e)Z{f), where / C e is drawn from the distribution 
P{f\e). However, this is true only if the coset e + A is cleanable. More generally, the initial state may 
have an error X{e)Z{g). Since Z-type errors commute with Tail, the hnal error becomes X{e)Z{g+f), 
where / C e is a random vector as above. 

The above discussion implicitly assumes that the initial error X{e) already includes a recovery 
operator applied by the decoder to correct the preceding memory errors. Ideally, the recovery cancels 
the error modulo stabilizers. In this case X (e) itself is a stabilizer and the transversal gate Taii creates 
no additional errors. If X{e) is not a stabilizer but the coset e + A is cleanable, application of Taii can 
possibly create new errors Z{f) with f + e. However, all these new Z-errors are either detectable 
or trivial due to the cleanability assumption. Thus, if the decoder has a sufficiently small list of 
candidate initial errors X{e), it might be possible to correct the new Z-errors at the later stage. On 
the other hand, if the coset e + A is not cleanable, at least one of the new Z-errors is undetectable 
and non-trivial, so the the encoded information is lost. Although the decoder cannot test cleanability 
condition, this can be easily done in numerical simulations where the actual error is known at each 
step. A good strategy is to test the cleanability condition prior to each application of the logical 
T-gate and abort the simulation whenever the test fails. Conditioned on passing the cleanability 
test at each step, the only effect of the transversal gates Tau is to propagate pre-existing X-errors to 
new Z-errors as described above. In this way the simulation can be performed using the standard 
stabilizer formalism, even though the circuit may contain non-Clifford gates. 

6 Maximum likelihood error correction and gauge fixing 

In this section we develop error correction and gauge hxing algorithms for simulation of logical 
Chfford-|-T circuits in a presence of noise. A fault-tolerant protocol implementing such logical circuit 
is modeled as a sequence of the following elementary steps: 

• Memory errors, 

• Noisy syndrome measurements, 

• Code deformations extending the gauge group to a larger group, 

• Code deformations restricting the gauge group to a subgroup, 

• Transversal logical Clifford gates, 

• Transversal logical T-gates. 


15 



We shall label steps of the protocol by an integer time variable t = 0 , 1 ,... ,L. At each step t the 
logical state is encoded by a quantum code of CSS type with a stabilizer group St and a gauge group 
Qt dehned as 

St = CSS{At,Bt) and Qt = CSS{Bt, A). 

Each code has one logical qubit with logical Pauli operators = W(l), = Z{1) and n physical 

qubits. A state of the physical qubits at the t-th step will be represented as EtUtEj, where ujt is 
some logical state of the code St and Et is a Pauli error. To simplify notations, we shall represent 
n-qubit Pauli operators X{a)Z{b) by binary vectors e G such that e = (a, b). Here we ignore the 
overall phase. Multiplication in the Pauli group corresponds to addition in F 2 ”. Each Pauli error 
e G F 2 ” is contained in some coset of the gauge group Qt, namely, e + Qt- The decoder described 
below operates on the level of cosets without ever trying to distinguish errors from the same coset. 
This is justihed since the logical state Ut is invariant under gauge operators, see Section 3. We 
shall parameterize cosets of Qt by binary vectors as follows. Let At and Bt be some hxed generating 
matrices of subspaces + (1) and Bt + {1). Dehne a matrix 


Ct 



(23) 


Then the coset of a Pauli error e = (a, 6) is parameterized by a vector / = CtC = {Bta,Atb)- Note 
that the matrix Ct has size q x 2n, where q = 2 + dim (St). Accordingly, the gauge group Qt has 2^* 
cosets. Furthermore, if CSS{At,Bt) is a regular code, that is. At = Bt and Bt = At then At and Bt 
become generating matrices of B^ and A^ respectively. In this case Bta determines the coset of At 
that contains a while Atb determines the coset of Bt that contains b. Furthermore, q < n + 1 with 
the equality in the case of regular codes. 

We assume that memory errors occur at half-integer times t — 1 /2, where t = 1, .. ., L- A memory 
error is modeled by a random vector e G Fl"" drawn from some hxed probability distribution tt. For 
example, the depolarizing noise that independently applies one of the Pauli operators X, Y, Z to each 
qubit with probability p/3 can be described by 7 r(e) = nj^=i Cn+j), where 7 ri( 0 , 0) = 1 — p and 
7 ri(l, 0) = 7 ri( 0 ,1) = 7 ri(l, 1) = p/3. Let e* G F 2 ” be the memory error that occurred at time t — 1/2 
and -|- ... -f e* be the accumulated error at the t-th step. Assuming that e G Fl” is drawn from 
the distribution tt, the corresponding coset / = Cte is drawn from a distribution 


Ptif) = 

e: Cte=f 


(24) 


The coset of Qt that contains the accumulated error at the t-th step is 


/• = C,(e‘+e" + ... + e'). 


(25) 


Syndrome measurements occur at integer times t = 1 , 2 ,..., L. The measurement performed at 
the t-th step determines syndromes for some subset of bt stabilizers from the stabilizer group St- 
These stabilizers may or may not be independent. They may or may not generate the full group St- 
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The syndrome measurement reveals a partial information about the coset of the accumulated error 
/* which can be represented by a vector 


= Mtf G 

for some binary matrix Mt of size bt x q. We shall refer to the vector s* as a syndrome. A precise 
form of the matrix Mt is not important at this point. The only restriction that we impose is that 
MtCfC = 0 if e represents a logical Pauli operator, that is, e = (1,0) or e = (0,1). This ensures 
that syndrome measurements do not affect the logical qubit. In a presence of measurement errors 
the observed syndrome may differ from the actual one. We model a measurement error by a random 
vector e' E drawn from some hxed probability distribution Qt such that the syndrome observed 
at the t-th step is s* = Mtf^ + e'. The family of matrices Ct, Mt and the probability distributions 
Pt, Qt capture all features of the error model that matter for the decoder. The decoding problem can 
be stated as follows. 

ML Decoding: Given a list of observed syndromes ..., s^, determine which coset of the gauge 
group Ql is most likely to contain the accumulated error at the last time step. 

Below we describe an implementation of the ML decoder with the running time 0{Lc2^), where 
c = maxj Q. The decoder can be used in the online regime such that the syndrome is revealed 
only at the time step t. In this setting the running time is 0{c2Q per time step. The decoder also 
includes a preprocessing step that computes matrices Ct, Mt, the effective error model Pt,Qt, and 
certain additional data required for implementation of logical T-gates. The preprocessing step takes 
time in the worst case. Since this step does not depend on the measured syndromes, it can be 
performed offline. The preprocessing can be done more efficiently if the code has a special structure. 

Consider first the simplest case when all codes are the same and no logical gates are applied. Our 
model then describes a quantum memory with one logical qubit. Since all codes are the same, we 
temporarily suppress the subscript t in Qt, Ct, Mt, Pt, Qt and q, bt. For simplicity, we assume that 
the initial state at time f = 0 has no errors, that is, /° = 0. Consider some time step t and let / G F 2 
be a coset of Q. Let pt{f) be the probability that / contains the accumulated error at the t-th step, 
/ = C(e^ + ... + e*), conditioned on the observed syndromes ..., sL Since we assume that all 
errors occur independently, one has 

t t 

u=i u=i 

where we set /* = / and /° = 0. Also we ignore the overall normalization of Pt{f) which depends 
only on the observed syndromes. We set Po(/) = 1 if / = 0 and Po(/) = 0 otherwise. It will be 
convenient to represent the probability distribution pt{f) by a likelihood vector 

\Pi) = E '’<(/) l/> (27) 

/6Fi 
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that belongs to the Hilbert space of c qubits. Dehne c-qubit operators 


f,g&2 /eiFi 


(28) 


Then 

\pt) = iQtP)---iQ2P)iQiPm. 

Next we observe that P commutes with any X-type Pauli operator, 

(/|X(h)P|^ 7 ) = (/ + h\P\g) = P{f + g + h) = {f\P\g + h) = {f\PX{h)\g). 

Thus P is diagonal in the X-basis {H\f)}, where H is the Walsh-Hadamard transform on c qubits, 

H\!) = 2-‘PY.(-iro\g). 

9&2 


Note that P[^ = /. The above shows that a matrix P = HPH is diagonal in the Z-basis, 

f’l/) = P/)l/>. where P(f) = P(g). (29) 

Assuming that the likelihood vector pt-i has been already computed and stored in a classical memory 
as a real vector of size 2^, one can compute pt in four steps: 

1. Apply the Walsh-Hadamard transform H. 

2. Apply the diagonal matrix P. 

3. Apply the Walsh-Hadamard transform p[. 

4. Apply the diagonal matrix Qt. 

Obviously, steps (2,4) take time 0(2'^). Using the fast Walsh-Hadamard transform one can implement 
steps (1, 3) in time 0{c2^). Overall, computing all likelihood vectors pi,..., pl requires time 0{Lc2^). 
Finally, the most likely coset of Q at the hnal time step is = argmaxj piif)- It can be found in 
time 0(2'^). If one additionally assumes that the hnal syndrome measurement is noiseless, there are 
only four cosets consistent with the hnal syndrome s^. In this case the most likely coset can be 
computed in time 0(1). 

Less formally, the ML decoding can be viewed as a competition between two opposing forces. 
Memory errors described by P-matrices tend to spread the support of the likelihood vector over 
many cosets, whereas syndrome measurements described by Q-matrices tend to localize the likelihood 
vector on cosets / whose syndrome Mf is sufficiently close to the observed one. Typically, the most 
likely coset coincides with the coset of the accumulated error if the likelihood vector pt remains 
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sufficiently peaked at all time steps. Once the likelihood vector becomes flat, a successful decoding 
might not be possible and the encoded information is lost. 

Next let us extend the ML decoder to code deformations that can change the gauge group Qt, 
see Section 3. For concreteness, assume that a code deformation mapping Qt_i to Qt occurs at time 
t — 1/4. We assume that for any time step t one has either Qt-i ^ Qt or Qt-i ^ Qt- The general case 
can be reduced to one of these cases by adding a dummy time step with a gauge group Qt-i H Qti 
such that first Qt-i is reduced to Qt-i H Qt and next Qt-i H Qt is extended to Qt- 
Case 1: Qt-i Q Qt- Then each coset of Qt-i uniquely determines a coset of Qt- In other words, 
/* = Dtf*~^ for some binary matrix Dt of size q x Ct-i- Note that q < Ct-i since a larger group has 
fewer cosets. Then the code deformation is equivalent to updating the likelihood vector according to 

Ia>^ E ft(/)IA/>. (30) 

Let us show that the updated vector can be computed in time Indeed, for any fixed g G F 2 * 

a linear system Dtf = g has solutions / G One can create a list of all solutions / and 

sum up the amplitudes pt{f) in time This would give a single amplitude of the updated 

likelihood vector. Since the number of such amplitudes is 2'^% the overall time scales as 0(2^*“^). A 
good strategy for numerical simulations is to choose matrices Ct in Eq. (23) such that Ct is obtained 
from Ct-i by removing a certain subset of Ct-i — Ct rows. Then /* is obtained from by erasing 
a certain subset of bits and the update Eq. (30) amounts to taking a partial trace. 

Case 2: Qt-i ^ Qt- Note that q > q_i since a smaller group has more cosets. Then each coset of 
Qt-i splits into cosets of Qt- Thus = Dtf^ for some binary matrix Dt of size Ct-i x q. As 

we argued in Section 3, the coset /* is drawn from the uniform distribution on the set of all cosets of 
Qt contained in see Eq. (9). Then the code deformation is equivalent to updating the likelihood 
vector according to 

Ip,) ^ E a(A/)I/). (31) 

The same arguments as above show that the update can be performed in time 0(2‘^‘). A good 
strategy for numerical simulations is to choose matrices Ct in Eq. (23) such that Ct is obtained from 
Ct-i by adding a certain subset of q — Ct-i rows. Then /* is obtained from by inserting random 
bits on a subset of coordinates and the action of Dt on the likelihood vector amounts to taking a 
tensor product with the maximally mixed state of q — Ct-i qubits. 

Next let us discuss transversal logical gates. For concreteness, assume that a logical gate applied 
at a step t occurs at time t + 1/4, where f = 0, ...,L — 1. A logical gate can be applied at a step 
t only if the corresponding code is regular, that is, St = Qt = CSS{At, Bt), where Bt = At- We 
only allow logical gates from the Clifford+T basis and assume that the code St obeys transversality 
conditions of Lemma 1. Let V be one of the transversal operators Haii, Sail, Tail defined in Section 4. 

We start from logical Clifford gates. Let Pt be the Pauli operator that represents the accumulated 
error + ... + eL The initial state before application of V has a form PtOJtP} , where ujt is a logical 
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state of the code St. The hnal state after application of V is 

VPtUtP}V^ = QtCotQl 

where Ot = VutV^ is a new logical state and Qt = VPtV^ is a new Pauli error. Thus the likelihood of 
a coset PQt before application of V must be the same as the likelihood of a coset QQt after application 
of V, where P is an arbitrary Pauli error and Q = VPV^. Note that the coset QQt depends only 
on the coset of P since VQtV^ = Qt, see the proof of Lemma 1. Thus the action of V on cosets of 
Qt dehned above can be described by a linear map v : —)■ F^L For concreteness, assume that 

cosets / G F 2 are parameterized as in Eq. (23) such that an error X{a)Z{b) belongs to a coset 
/ = (a,/9), where a = Bta and Q = Ath. If V implements a logical H-gate then v{a,Q) = {Q,a). 
If V implements a logical F-gate then v{a,(3) = {a,a + (3). One can compute the action of any 
other Clifford gate on cosets in a similar fashion. This shows that application of V is equivalent to 
updating the likelihood vector according to 

ipt) ^ Pt(f)ivf}- (32) 


This update can be performed in time 0(2'^‘). 

Next assume that V = Tail is a transversal T-gate. As we argued in Section 5, it is desirable to 
correct pre-existing memory errors of A-type before applying the T-gate. The ML decoder performs 
error correction by examining the current likelihood vector pt and choosing a recovery operator from 
the most likely coset of X-errors 


a* = argmaxN pt{a,l3). (33) 

a. ' ^ 

/3 

Applying a recovery operator X{a*) is equivalent to updating the likelihood vector according to 

\pt) pt{a, (3)\a + a*, (3). 

a,p 

This update can be performed in time 0(2'^‘). 

Recall that for regular CSS codes a Pauli error (a, 6) belongs to a coset {a, (3) = Ct{a,b), where 
a = Bttt is the coset of At that contains a and [3 = Atb is the coset of Bt that contains b. Let Aft be a 
set of cleanable cosets of At, see Dehnition 3. As we argued in Section 5, a lookup table of cleanable 
cosets can be computed in time 0(2"') at the preprocessing step. The next step of the ML decoder 
is to project the likelihood vector on the subspace spanned by cleanable cosets, 

\pt) 'G- EE Pt{a,/3)\a,/3). 

aeAft 13 

This update can be performed in time 0(2'^‘). 
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Now the decoder is ready to apply the transversal gate Tail. Here we assume for simplicity that 
Tall = r®". For each cleanable coset a let e(a) G Fg be some hxed representative of a that satishes 
condition of Dehnition 3, that is, Bte{a) = a and the support of e{a) contains no vectors from 
A:^- n (9. A lookup table of vectors e{a) must be computed at the preprocessing step. As was argued 
in Section 5, a composition of Tail and a twirling map maps a pre-existing Pauli error X{e)Z{g) 
to a probabilistic mixture of Pauli errors X{e)Z{f -|- g), where / C e is a random vector drawn from 
distribution P(/|e) dehned in Eqs. (21,22). Since the twirling map has no effect on the likelihood 
vector, it can be skipped. Thus Tau updates the likelihood vector according to 

ip.>-E E + Atf). (34) 

“./3 fQe{a) 

Dehne an auxiliary operator Pq, such that 

r„l/3>= + (35) 

/Ce(«) 


This operator acts on the Hilbert space of m qubits, where m = 
Eq. (34) can be written as 


\Pt) ^ 



\Pt)- 


We note that Pa is diagonal in the X-basis since 


dim (H^) 


l-l-dim(^t). Then 


(36) 


r„= 5^ P(f\e{a))X(Af). 

fCe(a) 


Let H be the Walsh-Hadamard transform on m qubits. Dehne an operator Pa = HTaH which is 
diagonal in the Z-basis of m qubits. The diagonal matrix elements of Pa are 

(,3|f„|f)) = ^ P(f\e(a))(-irx/^ (37) 

/Ce(«) 

Below we consider some hxed a and use a shorthand e = e{a). Substituting the explicit expression 
for P(/|e) from Eq. (21) into Eq. (37) yields 

(/3|f„|/3)= 2-l"l ^(-1)H9+M/2+/3^^T. (38) 

geBt(e)nBt(e)^ /Ce 


The sum over / vanishes unless the restriction of the vector AJ(] onto e coincides with g. Dehne a 
diagonal n x n matrix that has non-zero entries in the support of e, that is, diag(Je) = e. Then 
we get a constraint g = JeAj(3 and g G Bt{e) fl This shows that 


(/S|fa|/?) 


(_l)iPeAT/3| if J^A]^CBtie)nBtie)^, 
0 otherwise. 


(39) 
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The action of the diagonal operator Tq, on any state of m qnbits can be computed in time 0(2™') 
using Eq. (39). Using the fast Walsh-Hadamard transform on m qubits one can compute the action 
of the original operator T^^ dehned in Eq. (35) on any state of m qubits in time 0(m2™). The number 
of cleanable cosets a is upper bounded by the total number of cosets of At- The latter is equal to 2^, 
where k = dim {Aj~) = 1 + dim {Bt). The updated likelihood vector can be computed by evaluating 
each term a in Eq. (36) individually and summing up these terms. This calculation requires time 
= 0(ct2^‘). Overall, application of the transversal T-gate requires time 0{ct2^^). 

7 Logical Clifford+T circuits with the 15-qubit code 

In this section we apply the ML decoder to a specihc fault-tolerant protocol which similar to the 
one proposed by Anderson et al [14]. The protocol simulates a logical quantum circuit composed of 
Clifford and T gates by alternating between two error correcting codes: the 7-qubit color code (a.k.a. 
the Steane code) and the 15-qubit Reed-Muller code. We shall label these two codes C and T since 
they provide transversal Clifford gates and the T-gate respectively. Both C and T are regular CSS 
codes that can be obtained by the gauge hxing method from a single subsystem CSS code that we call 
the base code. This section may also serve as simple example illustrating the general construction of 
doubled color codes described in Sections 8-12. 

Our starting point is the 7-qubit color code. Consider a graph A = (U, E) with a set of vertices 
V = [7] and a set of edges E shown on Fig. 1(a). For consistency with the subsequent sections, we 
shall refer to the graph A as a lattice. Vertices of A will be referred to as sites. We place one qubit 
at each site of the lattice. Let f^, f"^, P T F 2 be the three faces of A, see Fig. 1(b), and 5 C F 2 be 
the three-dimensional subspace spanned by the faces, S = {f^, P, P)- The faces f^,P,P can also 
be viewed as rows of a parity check matrix 




1 

f 
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The 7-qubit color code has a stabilizer group CSS(5, S). This is a regular CSS code with six stabilizers 
A(/*), Z{p). The code distance is d = d{S) = 3. Minimum weight logical operators are Xp), Zp), 
where 

u = l + p = e^ + e^ + e^ 

Recall that e* denotes the standard basis vector of the binary space. The subspace S is doubly even, 
1/1=0 (mod 4) for all / G 5. Lemma 1 implies that the color code has transversal logical gates ff 
and S that can be realized by operators Haii = and Sail = 5'®'^. 

Suppose now that each site of the lattice A contains two qubits labeled A and B. Let us add one 
additional qubit labeled C. We assume that C is placed next to the lattice A such that C and u 
form one additional face, see Fig. 1(c). This dehnes a system of 15 qubits that are partitioned into 
three consecutive blocks, [15] = ABC, where |A| = \B\ = 7 and ICI = 1. 
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(a) 


(b) 


(c) 


Figure 1: (a) Color code lattice A. The 7-qubit color code is defined by placing one qubit at each 
site of A. (b) Faces /^/^,/^ and the subset u define stabilizers X(/*), Z’(/*) and logical operators 
X{u), Z{u) respectively, (c) A doubled color code is obtained by making two copies of A labeled A 
and B placed atop of each other, and adding one additional qubit labeled C. 


The T-code that provides a transversal T-gate has a stabilizer group CSS(T, T), where T C 
is a four-dimensional subspace spanned by double faces of the lattice A and by the all-ones vector 
supported on the region BC, 

T={r{A] + f‘lB]} + {BC}. (40) 

Here i = 1,2,3 and BC = 1[BC]. One can also view basis vectors of T as rows of a parity check 
matrix 



Here the first three rows stand for /*[A] -|- f^[B] and the last row stands for BC. A direct inspection 
shows that T is triply even, |/| = 0 (mod 8) for all / G T. By Lemma 1, the T-code has a 
transversal logical T-gate realized by an operator Taii = T®^®. Let us explicitly describe Z-stabilizers 
of the T-code. Define a subspace Q C spanned by double edges of A, 

g = {l[A] + l[B] : leE). (41) 

For example, a double edge I = (2,3) E E gives rise to a basis vector l[A] +l[B] = -|- -|- e® -|- 

see Fig. 1(c). Clearly, double edges l[A] + l[B] have an even overlap with double faces P[A] + P[B] 
as well as with the vector BC. This shows that Q <Z T. Define also a subspace C C F 2 ® spanned by 
single faces of A, including the extra face formed by u and C, namely 

C = (/‘M) + (/‘lB|) + {^lB|+C), (42) 

where i = 1, 2, 3 and C = 1[C]. A direct inspection shows that C CT and, moreover, 

f = C + g. (43) 
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To summarize, Z-stabilizers of the T-code fall into three classes: (i) edge-type stabilizers Z{f), where 
/ = /[^] + l[B] is a double edge of the color code lattice A, (ii) face-type stabilizers Z{g), where 
g = /*[A] or = /*[-B] is a face of the lattice A, and (hi) a special face-type stabilizer g = u[B] C 
that represents the extra face connecting C and u. These stabilizers are not independent. For 
example, the product of edge-type stabilizers over any closed loop on the color code lattice A gives 
the identity. Minimum weight logical operators of the T-code can be chosen as X{A) and Z{u:[A\). 
The code distance is d = 3 since d{T) = 3 and d{T) = 7. 

The (F-code that provides transversal Clifford gates has a stabilizer group CSS(C,C), where C 
is dehned by Eq. (42). Stabilizers of this code can be partitioned into three classes: (i) stabilizers 
X{P[A]), Z{f^[A]) dehne the 7-qubit color code on the region A with logical operators X{uj[A\), 
Z{uj[A\), (ii) stabilizers X(/*[i?]), Z(/®[T]) dehne the 7-qubit color code on the region B with logical 
operators X{uj[B]), Z{u:[B]), and (c) stabilizers X{u:[B] -|-C), Z{u:[B] -\-C) dehne the two-qubit EPR 
state |0, 0) -|- |1,1) shared between B and C such that the hrst qubit of the EPR state is encoded by 
the 7-qubit color code. Thus any logical state of the C-code has a form 

IV’i) = («|0 l) + P\Il))a ® (|0 l0) + |1l1))bc, (44) 

where a,/3 G C are some coefficients and |0 l), |1l) are the logical basis states of the 7-qubit color 
code. By discarding qubits of BC one can convert the C-code into the 7-qubit color code. Thus these 
two codes have the same logical operators, the same distance, and the same transversality properties. 
In particular, the C-code provides a transversal implementation of the full Clihord group. This can 
also be seen from Lemma 1 by noting that C = C and C is doubly even, |/| = 0 (mod 4) for all f &C. 

The base code appears as an intermediate step in the conversion between C and T codes. The 
stabilizer group of the base code is dehned as the intersection of stabilizer groups CSS(C,C) and 
CSS(T, T) describing the C and T codes. Thus the base code has a stabilizer group CSS(C fl T, C fl T). 
By dehnition, C C T, see Eq. (43), that is, C fl T = C. We claim that T ^C. Indeed, all double face 
generators of T are contained in C by dehnition, so we just need to check that BC G C. Using the 
identity p + u = 1, see Fig. l(a,b), one gets 

BC = B + C= {p-^u)[B] +C = p[B] + {uj[B] + C) G C. (45) 

We conclude that the base code has a stabilizer group CSS(T, C) and the gauge group CSS(C, T) = 
CSS(C, T). The base code has distance d = 3 since d(T) = d{C) = 3. All three codes have the same 
logical operators dehned in Eq. (4). We summarize dehnitions of the three codes in Table 1. 

We are now ready to describe a fault-tolerant implementation of a logical circuit in the Clihord-I-T 
basis. Our protocol consists of an alternating sequence of rounds labeled C and T, see Fig. 2. Each C- 
round is responsible for measuring syndromes of all face-type stabilizers of the O-code, that is, X{f) 
and Z{f), where / is one of the vectors /*[A], /*[i?], or u}[B] -\- C. Each T-round is responsible for 
measuring syndromes of all edge-type stabilizers of the T-code, that is, stabilizers Z{l[A] -\-l[B]) with 
I G E. All measured syndromes are sent to the decoder. Typically, but not always, a logical Clihord 
gate (T-gate) is applied after each C-round (T-round). Whether or not a logical gate is applied 
depends on the outcome of a certain test that we call a syndrome test. A decoder is responsible for 
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Transversal gates 

X-stabilizers 

Z-stabilizers 

stabilizer group 

gauge group 

C-code 

Clifford group 

single faces 

single faces 

CSS(C,C) 

CSS(C,C) 

T-code 

T gate 

double faces 
BC 

single faces 
double edges 

css(r, f) 

css(r, f) 

Base code 


double faces 
BC 

single faces 

css(r,c) 

css(c,r) 


Table 1: The family of codes used in the protocol. 


choosing a recovery operator R which is applied at the end of every pair of C, T rounds passing the 
syndrome test. In the beginning of each round the decoder performs a code deformation such that 
the logical qubit is encoded by the C-code (T-code) in every C-round (T-round). Although the base 
code does not explicitly appear in the protocol, it is used by the decoder as an intermediate step 
in the code deformation, see Section 6. Namely, at the beginning of each C-round the gauge group 
changes according to 

css(r,r) ^ css(c,r) ^ css(c,c). 

At the beginning of each T-round the gauge group changes in the reverse direction. 

Let us now describe the syndrome test. Recall that the syndromes of X(/) and Z{f) are denoted 
^(/) and C(/) respectively, see Section 3. Consider some fixed pair of rounds C, T and let U be the 
logical Clifford gate applied in the C-round (set C = / if no logical gate have been applied). Let ^(/) 
and ({f) be the face-type syndromes measured in this round. Measuring the syndrome of a stabilizer 
Z{f) after application of U is equivalent to measuring the syndrome of P{f) before application of 
U, where P = UZW. Suppose P ~ X{a)Z{b), where a,b E F 2 . Dehne an updated syndrome 

Cc/(/) = <(/) + &((/)• 

Thus Cu{f) determines the syndrome of a stabilizer Z{f) that would be observed in the absence of 
the logical gate U. Let C(^[^] + ^[B]) be the edge-type syndromes measured in the T-round. We say 
that the pair of rounds C, T passes a syndrome test if 

C((|.41 + i|B]) + C(r[.4] + l‘\B]) + (uirw) + &(/‘|Bl) = 0 (46) 

for any face /* and for any pair of edges 1,1' E E such that I +1' = /*. In other words, I and I' are 
the two opposite edges forming the boundary of /*. 

In the absence of errors the syndrome test is always passed since the product of edge-type stabiliz¬ 
ers Z{l[A] + l[B]), Z{l'[A]+l'[B]) and face-type stabilizers Z(/*[A]), Z{f''[B]) equals the identity. Our 
protocol performs the syndrome test after each T-round, see Fig. 2. If the syndrome test fails, an ad¬ 
ditional pair of rounds C, T is requested and the process continues until the syndrome test is passed. 
We note that the syndrome test can fail for at least two reasons. First, any single-qubit W-error on 
some qubit j E AB that occurs inside the chosen T-round flips the syndromes C(^^] + ^[B]) on all 
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C-round 


T-round 


Optional rounds. 

Repeat until syndrome test is passed. 


Figure 2: A fault-tolerant implementation of a Clifford-|-T circuit. Red circles represent depolarizing 
memory errors and syndrome measurement errors. Transversal Clifford and T-gates are shown by C 
and T boxes. The boxes W and R represent the twirling map and the recovery operator respectively. 
Measurement boxes represent (partial) syndrome measurements. Face-type stabilizers of the C-code 
are measured in every C-round. Edge-type stabilizers of the T-code are measured in every T-round. 
Optional rounds are added only if the full syndrome of the T-code inferred from these measurements 
fails to pass a consistency test. The rounds continue in the periodic fashion. On average, the protocol 
applies one logical gate per round. 


edges I incident to the site u E A that contains j without flipping any face-type syndromes (because 
the latter have been measured before this error occurred). This would violate at least one constraint 
in Eq. (46). Secondly, any single measurement error for edge-type stabilizers Z{l[A] + l[B]) in the 
T-round or any single measurement error for face-type stabilizers that contribute to the updated 
syndromes Cu{f) would violate at least one constraint in Eq. (46). The purpose of the syndrome test 
is to ensure that neither of these possibilities occurs before asking the decoder to perform a recovery 
operation. 

Combining syndromes measured in any consecutive pair of rounds C, T provides the full syndrome 
for the T-code (in the absence of errors). Indeed, the syndrome of X{T) can be inferred from the 
syndromes of X{C) and Z{C) measured in the C-round since T T C and X{T) commutes with all 

operators measured in the T-round. Likewise, the syndrome of Z{T) can be obtained by combining 
the syndromes of X{C) and Z{C) measured in the C-round and the syndrome of Z{Q) measured in 

the T-round, since T = C + Q. By spreading the syndrome measurement for the T-code over two 
rounds we were able to keep the number of measurements per qubit in any single round reasonably 
small, which might be important for practical implementation. 

A transversal Clifford gate is applied after each C-round provided that the latest syndrome tests 
was successful. We performed simulations for a random Clifford-I-T circuit such that each Clifford 
gate is drawn from the uniform distribution on the Clifford group. 

A transversal T-gate is applied at the end of each T-round that passes the syndrome test. It 
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is preceded by a Pauli recovery operator R classically controlled by the decoder, see Section 6, and 
followed by a twirling map Wr that applies a randomly chosen stabilizer X{f) with f & T, 

^t(p) = T^,J2x{f)pX(f). ( 47 ) 

' ' /er 

Each round includes memory and syndrome measurement errors. We model memory errors by the 
depolarizing noise with some error rate p, that is, each qubit suffers from a Pauli error X, Y, Z with 
probability p/3 each. Within each round a memory error occurs before the syndrome measurement, 
see Fig. 2. We model a noisy syndrome measurement by an ideal measurement in which the outcome 
is flipped with a probability p. 

At any given time step the protocol can be terminated depending on the outcome of two tests: (1) 
logical error test and (2) cleanability test. A logical error test is performed at the end of each round 
by computing the most likely coset of errors consistent with the current syndrome. The test is passed 
if the most likely coset contains the actual memory error. A cleanability test is performed after each 
recovery operation. If E ~ X(a)Z(b) is the residual error left after the recovery, the test is passed if 
a + T is a cleanable coset, see Eq. (40) and Dehnition 3. We found that T has 996 cleanable cosets. 
The protocol terminates whenever one of the two tests fails. Accordingly, the number of logical gates 
implemented in the protocol is a random variable. Conditioned on passing the cleanability test, a 
transversal T-gate is implemented using the method of Section 5. The quantity we are interested in 
is a logical error rate dehned as pi = ^/g, where g is the average number of logical gates implemented 
before the protocol terminates. Here g includes both Clifford and T gates. 

The logical error rate pi was computed numerically using a simplihed version of the ML decoder 
that we call a sparse ML decoder (SMLD). It follows the algorithm described in Section 6 with two 
modihcations. First, SMLD models memory errors by a sparse distribution tt'^ approximating the 
exact distribution tt. We have chosen 7r^(e) = 7r(e) if e is a Pauli error acting non-trivially on at 
most one qubit and 7r^(e) = 0 otherwise (here we ignore the normalization). Replacing tt by vr^ in the 
update rules of Section 6 one can see that the matrix Pt dehned in Eqs. (24,28) becomes sparse and 
the update \pt) <(— Pt\pt) can be performed using sparse matrix-vector multiplication avoiding the 
Walsh-Hadamard transforms. The actual memory errors in the Monte Carlo simulation are drawn 
from the exact distribution tt. This approximation is justihed since none of the codes used in the 
protocol can correct memory errors of weight larger than one. Secondly, SMLD performs a truncation 
of the likelihood vector pt after each round in order to keep pt sufficiently sparse. The truncation was 
performed by normalizing pt and setting to zero all components with pt{f) < e, where e = 10“® is an 
empirically chosen cutoff value. We observed that for large error rates {p ~ 1%) the two versions of 
the decoder achieve the same logical error probability within statistical huctuations. On the other 
hand, SMLD provides at least lOx speedup compared with the exact version and enables simulation 
of circuits with more than 10,000 logical gates. Our results are presented on Fig. 3. We observed a 
scaling Pi = Cp^ with C ~ 182. Assuming that a physical Clifford-|-T circuit has an error probability 
p per gate, the logical circuit becomes more reliable than the physical one provided that pi < p, 
that is, p < Po = C~^ ~ 0.55%. This value can be viewed as an “error threshold” of the proposed 
protocol. Generating the data shown on Fig. 3 took approximately one day on a laptop computer. 
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Figure 3: Monte Carlo simulation of logical Clifford+T circuits with the sparse ML decoder (SMLD). 
Here p is the physical error rate and pl is the logical error rate dehned as pl = l/f/, where g is the 
average number of logical gates implemented before the protocol terminates. Each average was 
estimated using 400 Monte Carlo trials. 

8 Doubled color codes: main properties 

Our goal for the rest of the paper is to generalize the 15-qubit codes described in Section 7 to higher 
distance codes that can be embedded into a 2D lattice such that all syndrome measurements required 
for error correction and gauge hxing are spatially local. This section highlights main properties of the 
new codes. Let d = 2f -|- 1 be the desired code distance, where f > 1 is an integer. For each t we shall 
construct a pair of CSS codes labeled C and T that encode one logical qubit into rit = 2t^ + 8t^-|-6f —1 
physical qubits. These codes have transversal Clifford gates and the T-gate respectively. The codes 
are defined on the 2D honeycomb lattice with two qubits per site such that the gauge group of the 
C-code has spatially local generators supported on faces of the lattice. Most of these generators are 
analogous to face-type stabilizers in the 15-qubit example, see Section 7. There are also additional 
generators of weight two that couple pair of qubits located on the same face. The latter have no 
analogue in the 15-qubit example. The C-code can be converted to the regular color code on the 
honeycomb lattice by discarding a certain subset of qubits. Although the T-code does not have local 
gauge generators, it can be obtained from the C-code by a local gauge hxing. It requires syndrome 
measurements for weight-four stabilizers supported on edges of the lattice. The latter are analogous 
to edge-type stabilizers in the 15-qubit example. We shall also dehne a base code that appears as an 
intermediate step in the conversion between C and T codes. All three codes have distance d = 2t + l. 

The codes will be constructed from a pair of linear subspaces Ct,Tt F 2 * that satisfy 

TtCCtCCtCft (48) 
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and 


d{Tt) — 2t + 1. (49) 

The subspaces Ct and % have a special symmetry required for transversality of logical gates, namely, 
Ct is doubly even and % is triply even with respect to some subsets of qubits C [n*] such that 
\M^\ — \ M^\ = 1. Furthermore, Ct and 7* have spatially local generators (basis vectors) supported on 
faces of the lattice. The subspaces Ci and 7i coincide with C and T dehned in Section 7. Dehnitions 
of the three codes are summarized in Table 2. We shall refer to the family of codes dehned in Table 2 
as doubled color codes since each of them is constructed from two copies of the regular color code. 



Transversal gates 

Stabilizer group 

Gauge group 

C-code 

Clifford group 

CSS{Ct,Ct) 

CSS{Ct,Ct) 

T-code 

T gate 

css{Tt,rt) 

CSSiT,ft) 

Base code 


CSSiT,Ct) 

CSS{Ct,ft) 


Table 2: A family of 2D doubled color codes that achieves universality by the gauge hxing method. 
Each code has one logical qubit, n* = 2t^ + 8t‘^ + 6t — l physical qubits, and distance d = 2t + l. The 
subspaces Ct and Tt have spatially local generators supported on faces of the lattice. These subspaces 
are doubly even and triply even respectively. 

We expect that the new codes can be used in a protocol analogous to the one shown on Fig. 2 
to implement fault-tolerant Clifford-|-T circuits in the 2D architecture. The C-round of the protocol 
would be responsible for measuring gauge syndromes ^{f^) and Cif^) some complete set of 
generators Z^,...,/'" G Ct- A transversal Clifford gate can be applied after each C-round. The 
T-round would be responsible for measuring stabilizer syndromes some set of generators 

g^,... ,g’^ G Tt such that T = (/^,..., /”*, g^,..., g^). We will see that such generators g^^ can be 
chosen as double edges of the lattice, by analogy with the 15-qubit example. A transversal T-gate can 
be applied after each T-round. In the case of higher distance codes the C-round could be repeated 
several times to ensure that any combination of t syndrome measurement errors is correctable. We 
expect that the T-round requires less (if any) repetitions since the edge-type stabilizers measured 
in this round are highly redundant. An explicit construction of a fault-tolerant protocol that would 
suppress the error probability per gate from p to 0(p*^^) is an interesting open problem that we 
leave for a future work. 


9 Regular color codes 

The starting point for our construction is the standard color code on the hexagonal lattice [23]. For 
a suitable choice of boundary conditions such code has exactly one logical qubit and corrects t single¬ 
qubit errors, where t > 0 is any given integer (recall that the code distance is d = 2t -|- 1). Examples 
of the color code lattice for t < 3 are shown on Fig. 4. More formally, consider a lattice At such that 
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the lattice sites are triples of non-negative integers j = (ji, 72 , 73 ) satisfying 71 -I -72 -I -73 = 3t. Clearly, 
At is the regular triangular lattice. Note that 72 —71 = 73 —72 = 7 i —73 (mod 3) for any j G At. For 
each 6 G Z 3 consider a sublattice 

= { j G Aj : j 2 -ji = b (mod 3) }. (50) 

Dehne a color code lattice as 

Ai = A°UA?. (51) 

In other words, sites of A^ are triples of non-negative integers j = ( 71 , 72 , 73 ) satisfying ji -I- 72+73 = 3t 
and 72 —71 7 ^ 1 (mod 3). Sites of the triangular lattice At that are not present in A^ become centers 
of faces of A^. More formally, a subset / C At is called a face iff there exists j G A) such that / is 
the set of nearest neighbors of j in the triangular lattice At. A direct inspection shows that any face 
of At consists of four or six sites. Although the color code lattice is usually equipped with a face 
3-coloring, we shall not specify the colors since they are irrelevant for what follows. 




Figure 4: Color code lattices Ai, A 2 , and A 3 . Qubits are located at the lattice sites (empty circles). 
Each face / gives rise to a pair of stabilizer generators X{f) and Z{f). The color code dehned on 
the lattice A* has one logical qubit and corrects t single-qubit errors (code distance is d = 2 t -|- 1 ). 

The lattice A^ contains 

mt = \At\ = + 3t + 1 (52) 

sites and {rrit — l)/2 faces. We will use a convention that Aq is a single site. Let St Q F™* be 
a subspace spanned by all faces of At (recall that we identify a face / and a binary vector whose 
support is /). Note that St is self-orthogonal, since any pair of faces overlap on even number of sites. 
We shall need the following well-known properties of St, see Refs. [23, 18, 24]. 

Fact 1. The subspace St is doubly-even with respect to the subsets A^ and A^. Furthermore, 

|A?| - |A?| = 1. (53) 
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The orthogonal subspace is given by 


nS = St and no = St + 1, 


(54) 


Finally, 


d{St) = 2t + l. 


(55) 


For the sake of completeness, let us prove the hrst claim. 


Proof. We shall use shorthand notations S = St, m = rut, and = A^. Consider any stabilizer 
f E S. Then / is a linear combination of faces. One can always choose the numbering of faces such 
that / = Yli=i /* some k < s. We shall use an identity 




k 

5^ If/1-2 (mod 4) 

i=l 


(56) 


which holds for any vectors g^,... ,g^ G F™. Choosing S'* = /* O A^ one gets 


k 

|/nA"|-|/nA"| = ^(|/'nA»|-|/‘nA=|)-2 ^ (|/‘n/-'n A"| - l/'nf n A"|) (57) 

i=l 

modulo four. Consider any site j G A^. The nearest neighbors of j in the triangular lattice belong 
to either A° or A^ as shown on Fig. 5. By examining this hgure one can easily check that any edge 
of the color code lattice Kt connects some site of A*^ and some site A^. 


Q P 

A^ O OA° 

AO O 0^2 


Figure 5: A local neighborhood of a site j G A^. This shows that any edge of the color code lattice 
At = A° U A^ has one endpoint in A° and the other endpoint in A^. 

It follows that I/* n A°| — I/* n A^l = 3 — 3 = 0 for hexagonal faces /* and |/* fl A°| — |/* n A^j = 
2 — 2 = 0 for square faces /*. Likewise, consider any pair of faces /*, such that /* fl ^ 0. Then 
/* n p is an edge of the color code lattice At, so that 

If n f n A°| - If n f n a2| = 1 - 1 = 0. 

This proves that all terms in Eq. (57) are zero, that is, S is triply-even as promised. 

□ 
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Finally, let and oj^ be the vectors supported on the right, on the left, and on the bottom 

sides of the triangle that forms the boundary of A*. In other words, 

= : j- = 0}, * = 1,2,3. (58) 

Note that u] are logical operators, that is, ul G fl O. Furthermore, \u;l\ = 2t + 1 = d, that is, ojI 
are minimum weight logical operators. 


10 Doubling transformation 

Let us now describe a general construction of triply even subspaces inspired by Ref. [20]. Consider 
a pair of integers m, n and a pair of subspaces 


5 C F™ and T C F”. 

Let k = 2m + n. Partition the set of integers [k] into three consecutive blocks, [k] = ABC, such that 
|A| = \B\ = m and \C\ = n. Dehne a subspace W C F^ such that 

U = (/[A] + f[B] : / G 5) + T[C] + {BC). (59) 


Here we used a shorthand notation BC = 1[BC]. This dehnition can be rephrased in terms of the 
generating matrices of S and T. Recall that S' is a generating matrix of a linear subspace 5 if 5 is 
spanned by rows of S. Suppose S and T have generating matrices S and T respectively. Then U 
has a generating matrix 

\ S S 1 


U 


T 
T T 


(60) 


where the three groups of columns correspond to A,B, and C respectively. Here we only show the 
non-zero entries of U. The mapping from S and T to U will be referred to as a doubling map since 
it involves taking two copies of S. We shall use an informal notation U = 2S + T to indicate that U 
is obtained from S and T via the doubling map as described above. 

Given a subset M C [m] and a subset N (Z [n] let MMA^ C ABC be a subset obtained by 
choosing the subset M in the blocks A, B and choosing the subset N in the block C. 


Lemma 2. Assume S is doubly even with respect to some subsets C [m] and T is triply even 
with respect to some subsets C [n] such that 


\M+\-\M-\ + \N-\-\N+\ = 0 (mods). 


(61) 


Then U = 2S + T is triply even with respect to subsets . 


32 




Proof. Consider first a vector 


li = /M + flB] + g[C] e U, 

where f E S and g E T. By assumption, 

1/n M"*"! — 1/n M“| = 0 (mod 4) and fl fl iV“| = 0 (mod 8). (62) 

The identity \h fl K^ \ = 2\f fl M^| + l^? fl then implies 

Ihn/sT+l - |hn/s:-| = 0 (mods). 

Consider now a vector 

h' = h + BC eU. 

Then 

\h' n K^\ = 1/ n M^l + (|M±| - 1/ n M^l) + {\N^\ -{gn A^^l) = \M^\ + \N^\ -\gn iV^|. 
Taking into account Eqs. (61,62) one arrives at 

\h'r]K+\-\h'nK-\=0 (mods). 

Since any vector of U can be written as h or h', the lemma is proved. □ 

Next let us compute the orthogonal subspace U^. We specialize to the case when m,n are odd, 
whereas all vectors in S and T have even weight. This will be the case for applications considered 
below, where S and T dehne stabilizer groups of CSS-type quantum codes. 

Lemma 3. Suppose n and m are odd. Suppose S S and T ^ £. Then 

= if [A] + f[B] : f e£)+S^[A]+ f[C] + (BC). (63) 

Furthermore, 

U = if [A] + f[B] : fE£)+ S[B] + f[C] + (BC). (64) 

and 

dipt) = min {d{S), d{T) + 2}. (65) 

Proof. Consider an arbitrary vector h By dehnition, the inclusion h G is equivalent to 

hA + hsES^, hcET^, and hB®hcE£. (66) 

A direct inspection shows that Eq. (66) holds if h belongs to each individual term in Eq. (63) which 
proves the inclusion 3 in Eq. (63). Conversely, suppose h E lA^. We have to prove that h is contained 
in the sum of the four terms in Eq. (63). Set f = hs A hs E £ and f = hs + i A hs E O. In both 
cases f E £. Replacing h hj h + f[A] + f[B] we can make hs = 0 or hs = 1. Since 1 G 5-*“, the hrst 
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condition in Eq. (66) implies Ha € iS-*-. If hs = 1, replace h hj h + 1[B] + 1[C]. This makes hs = 0 
and does not change the second condition in Eq. (66) since 1 G T^. The last condition in Eq. (66) 
then implies he E S. We conclude that h = g[A] + g'[C] for some g G and g' G T. This proves 
Eq. (63). 

To prove Eq. (64) we note that odd-weight vectors in can only originate from the term iS'*' [A\. 
Restriction to the even subspace replaces this term by 5[A]. Since we already know that U contains 
all vectors f[A\ -|- f[B] with / G T, we can move S from A io B without changing U. This proves 
Eq. (64). 

It remains to prove Eq. (65). Choose any vectors f* G fl O and g* G fl O such that 
d{S) = \f*\ and d(T) = \g*\. Choose any i G [m] and consider vectors x = e''[A\ + e*[i?] -|- g*[C] and 
y = f*[A]. A direct inspection shows that x^y & fl O. Therefore d{U) < |a:| =2+\g*\ = d{T) -|- 2 
and d(U) < \y\ = \f*\ = d{S). This proves the inequality < in Eq. (65). Let us prove the reverse 
inequality. Consider any vector h G fl O. It must satisfy Eq. (66). Since h G O, the condition 
hs ® he G £ implies G O. Consider two cases. Case 1: he & O. Then the second condition 
in Eq. (66) implies he G fl (9, that is, \he\ > d{T). Furthermore, since h,he,hA G we infer 

that hB & O which implies \h\ = \hA\ + Ih^l -|- \he\ > 2 + d{T). Case 2: he G £. Then the 

condition hs ® he ^ £ implies hs G T, that is, hA + hs & O. The hrst condition in Eq. (66) implies 
hA + hs G n C>, that is, > d{S). By triangle inequality, |h^| -|- |hs| > d{S) and thus 

\h\ > d{S). This proves the inequality > in Eq. (65). □ 

The above lemma has the following obvious corollary. 

Corollary 1. Suppose n andm are odd, S ®£, andT is self-orthogonal. ThenU is self-orthogonal. 

Proof. Self-orthogonality of T implies B P T. Thus the hrst, the second, and the third terms of 
Eq. (59) are contained in the hrst, the third, and the fourth terms of Eq. (63) respectively. □ 


11 Doubled color codes: construction 

Let St be the subspace spanned by faces of the color code lattice constructed in Section 9. Recall 
that St ® F™* where rut = [A^j = 3f^ -|- 3t -|- 1. A doubled color code with distance d = 2t -\-l will 
require Ut physical qubits, where uq = 1 and 

nt = 2mt + nt-i for f > 1. (67) 

Solving the recurrence relation gives 

nt = 2t^ + 6t^ -|- 6f -|- 1. (68) 

For instance, rii = 15, 77-2 = 53, and = 127. Dehne a family of subspaces Tt ® F 2 * such that 
To = (0) C F 2 and 

Tt = 2St + Tt-i, (69) 
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In other words, % is obtained by applying the donbling map of Section 10 with S = St and T = Tt-i- 
To describe this more explicitly, partition the set of integers [rzi] into 2t + 1 consecntive blocks as 

[nt\ = AtBt... ^ 2 - 82 ^ 1 - 81^0 where Ar = Br = A^. (70) 


In other words, Ar and Br represent two copies of the color code lattice A^. Then 


t t 

Tt = Y.{f[Ar] + f[Br\ fe Sr) + 

r=l r=l 


(71) 


We can also describe % by its generating matrix. Snppose St is a generating matrix of St snch that 
rows of St correspond to faces of the color code lattice Kf Then % has a generating matrix 


Ai Bt 

At-i Bt-i 


A2 -82 

Ai Bi 

Aq 

St St 

T 

St-i St-i 

T 

T 

1 

T 

82 82 

T 

T 

81 81 

T 

T 

1 


Here the hrst line indicates which block of qubits contains a given group of columns. We shall use 
the subspace Tt to construct the T-code as dehned in Table 2 . Let us prove that % has the properties 
stated in Section 8 . 


Lemma 4. The subspace % is triply even with respect to some subsets C [nt] satisfying 


liv,-*'! - lA'.l = 1- 


(72) 


Proof. We shall use induction in t. The base of induction, t = 0, corresponds to uq = 1 and 
To = ( 0 ) C F 2 . Clearly, To is triply even with respect to subsets 

K = {1], N^ = (ll (73) 

which obey Eq. (72). Consider now an arbitrary t. We already know that St is doubly even with 
respect to the subsets and Aj such that 

|A?| - |A?| = 1 . (74) 
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see Fact 1. Define 


At=Bt = ^'t and A; = B; = Al (75) 

Here we consider Af and as subsets of At and Bt respectively. Choose 

Aff = AfBfNlt, t > 1. (76) 

Combining Eqs. (74,75,76) and assnming that |iVjti| — |-/V^'li| = 1 one gets 

|Af+| - lAfrl = 2(|A?| - |A?|) + - \Nlt\ = 1. (77) 

This proves Eq. (72) for all t > 0. Fnrthermore, Eqs. (72,74,75) imply 

i7i+i-|7i,-|-|A'tl + |A'rl = ') 

for all t > 0. This shows that condition Eq. (61) of Lemma 2 is satished for = Af = Bf. The 
lemma implies that Tt is triply even with respect to the snbsets for all t > 0. □ 

Let us use induction in t to show that % has distance 

d{Tt) = 2t + l. (78) 

Indeed, d{%) = 1 since = (O)-*- = F 2 and the only odd-weight vector in F 2 is 1. Furthermore, the 
color code on the lattice has distance 2t -|- 1, that is, d{St) = 2t + 1, see Fact 1. Assuming that 
d{Tt-i) = 2t —1 and using Eq. (65) of Lemma 3 we infer that d{Tt) = min {<7(5*), 2 -|- d{Tt-\)} = 2t-|-l 
which proves Eq. (78) for all t > 0. 

To construct the T-code we shall also need a subspace 7t, see Table 2. It will be convenient to 
rewrite the partition in Eq. (70) as [n*] = A^BtCt, where 

Ct = At-iBt-i ... AiBiAq. 

Applying Eq. (64) of Lemma 3 and taking into account that St = St, see Fact 1, one gets To = 0 and 

= (/[A] + f[Bt] ■ f & £) + St[Bt] + it-ilCt] + (BtAt-i) (79) 

for f > 1. Let u)t G S^ AO he some hxed minimum weight logical operator of the regular color code 
such that \ujt\ = 2t + l. Later on we shall choose Ut as dehned in Eq. (58). Then 1 = ojt + g ioi some 
g G St, see Eq. (54) of Fact 1. This implies BtAt-i = u)t[Bt\ + g[Bt\ + At-i. Since g[Bt\ is contained 
in the second term in Eq. (79), we can replace the last term by {uJt[Bt] -|- At-i). Likewise, 

At_i = ijJt-i[At-i] A f[Ct\ 

for some / G Tt-i- Here we noted that both l[At_i] and uJt-i[At_\\ are contained in T{!ii fl O, so 
that the sum of them is contained in Tt^i A £ = 7t_i. Since /[Cf] is contained in the third term in 
Eq. (79), we can rewrite the last term as {ut[Bt\ + uJt-i[At-i]). Therefore 

'Bt = (/[A] + f[Bt\ '■ / G F) -|- St[Bt] + 7t-i[Ct] + {ut[Bt] + Ut-i[At-i]) (80) 
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for t > 1. Note that here we can add a term iSj [A*] since f[At] + f[Bt] is contained in the hrst term in 
Eq. (80) for any f E St- We can also describe Tt by its generating matrix. Suppose Et is a generating 
matrix of the even subspace such that rows of Et correspond to edges of the color code lattice A^, 
that is, each row of Et has a form e“ + e'^ for some edge {u, v) of A^. Then Tt has a generating matrix 



We shall refer to the hrst two groups of rows as face-type generators and edge-type generators. 

To construct the C-code we shall need a subspace Ct C dehned as 

t t 

Ct = '^^{f[Ar] -\- g[Br] : f , g E Sr) ''^^{BrAr-l) . (81) 

r=l r=l 

By comparing Eqs. (71,79,81) one can see that 

Tt<TCt = Ct<Tft. (82) 

A generating matrix of Ct can be obtained from Tt by removing all edge-type generators. A direct 
inspection shows that generators of Ct can be partitioned into mutually disjoint subsets supported 
on regions 


Mt — At, Mt-i 


BtAt-i, Mt-2 — Bt-iAt-2, 


Mo — BiAo- 
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Thus we can analyze properties of Ct on each region Mr separately. Generators of Ct supported on 
Mt have a form /[Aj] with / G St- These generators describe the regular color code St. Generators 
supported on a region Mr describe the two-qubit EPR state |0, 0) -|- |1,1) shared between and 
Ar such that the two qubits are encoded by the regular color code 5^+1 and Sr- As a consequence 
we obtain 

Corollary 2. The subspace Ct is doubly even with respect to the subsets C At, see Fact 1. 
Furthermore, d{Ct) = 2t -|- 1. 

At this point we have proved all properties of the subspaces Ct, % stated in Section 8 except for 
the spatial locality. Let A be the honeycomb lattice with two qubits per site. We shall allocate a 
triangular-shaped region of A isomorphic to the color code lattice A^ to accommodate the blocks of 
qubits At and Bt- The two blocks can share the same set of sites since each site contains two qubits. 
By a slight abuse of terminology, we shall identify At and the region of A accommodating At and Bt. 
The regions A* and Ai_i are placed next to each other as shown on Fig. 6. 

Gonsider the generating matrix Tt dehned above. We shall say that a vector / G is spatially 
local if its support is contained in a single face of the lattice. By dehnition, face-type and edge-type 
generators of Tt are spatially local. Gonsider some row of Tt in the bottom group. It has a form 

UJr,r-l = UJr[Br] -|- UJr-l[Ar-]\, r = 1, . . . ,t. 

The above arguments show that we are free to choose different logical operators Ur in different 
generators Let us use this freedom to choose 

^r,r-l = ^\.[Br\ + Oji_i[Ar-\\ (83) 

where cj* is the minimum weight logical operator supported on the Ath boundary of the lattice A^, see 
Eq. (58). Furthermore, we can choose i and j such that the Ath boundary of A^ is located next to the 
j-th boundary of Ar_i, see Fig. 6. Then the generator ujr,r-i has a shape of a loop that encloses the 
free space separating the regions A^ and Aj.-i. Since |a;*| = 2r-|-l, we get \ujr,r-i \ = 2r-|-l-|-2r —1 = 4r. 
We shall explain how to reduce the weight of the generators uJr^r-i and make then spatially local in 
the next section. 

The above discussion also shows that a code deformation transforming the G-code to the T-code 
requires only spatially local syndrome measurements. Indeed, suppose pl is some logical state of the 
G-code, see Table 2. One can hrst apply a reverse gauge hxing that extends the gauge group of 
from GSS(Cf,Ci) to GSS(Ct,7t). This can be achieved by applying a random element of the group 
Z{Tt). This also restricts the stabilizer group of p^ from GSS(G,Ct) to GSS(7t, G). We note that 

ft = Ct+ Y. {l[A]+l[B]), 

l={u,v) 

where the sum runs over all edges of the sub-lattices Ai,..., A^. Thus a gauge hxing that extends 
the stabilizer group of pi from CSS{Tt,Ct) to GSS(7i, %) can be realized by measuring syndromes of 
the edge-type stabilizers Z{l[A] + l[B]) and applying a suitable recovery operator. 
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Figure 6 : Embedding of the doubled color code Ts into the hexagonal lattice A with two qubits per 
site. The lattice is divided into four disjoint regions Aq, Ai, A 2 , A 3 where At accommodates two copies 
of the color code St labeled At and Bt (not shown). Recall that Aq is a single site. To enable a more 
regular embedding we twisted one corner of each lattice At compared with Fig. 4. Sites represented 
by solid circles can be ignored at this point. To obtain the next code 71 one should attach A 4 to the 
bottom side of A3. To obtain Ts one should attach A5 to the left side of A4. The process continues in 
the alternating fashion. Gauge generators live on edges and faces of the lattice. There are also non¬ 
local gauge generators ujt,t-i of weight 4f connecting R-qubits on the boundary of A* and A-qubits 
on the boundary of At_i (not shown). 

12 Weight reduction 

In this section we transform the G-code into a local form. This requires two steps. First, we show 
how to represent each non-local gauge generator Ur^r-i, see Eq. (83), as a sum of spatially local 
generators of weight at most six and a single non-local generator of weight two. Secondly, we show 
how to represent each of the remaining non-local generators as a sum of spatially local generators 
of weight two. Each of these steps extends the code by adding several ancillary qubits and gauge 
generators. We add the same ancillary qubits and gauge generators to both C and T codes to preserve 
the local mapping between them. Accordingly, we have to prove that the extended versions of the 
C and T codes have the same distance and the same transversality properties as the original codes. 

We begin by setting up some notations. Below we consider some hxed value of t and r = 1,..., f. 
Let us denote the sites of A^ lying on the boundary facing A^-i as u^,u^,... The ordering 

is chosen such that is the “twisted” corner of A^, see Fig. 7. Let us denote the sites of A^-i 

lying on the boundary facing A^ as ..., The ordering is chosen such that u* and u* are 
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next-to-nearest neighbors, see Fig. 7. Using these notations one can rewrite Eq. (83) as 

2r+l 2r—l 

(^r,r-l — ^ ^ ^ ^ 'U^[A^_l]. (84) 

i=l i=l 

Consider now sites of A lying in the free space separating and A^-i. These sites are indicated by 

hlled circles on Figs. 6,7. Denote these sites as w^, w ^,..., see Fig. 7, and let 

Dr = . . . , 

The ordering is chosen such that tc® is a nearest neighbor of m® and f® for all 1 < i < 2r — 1. 
Furthermore, tc^® is a nearest neighbor of 


1 / 1/6 ( 


1/5 








v> 


1 / 1/5 • • 1 / 1/4 1 / 1/3 0 0 1 / 1/2 




Figure 7: Sites m® on the boundary of A 3 (blue), sites n® on the boundary of A 2 (red), and sites w® in 
the region (hlled circles). The lattice is rotated 60° counter-clockwise compared with Fig. 6 . 

Let us add 2r ancillary qubits such that each site of Dr contains one qubit. The total number of 
physical qubits in the lattice A becomes 

t 

Nt = nt + J2^^ = 2^^ + + 7t + 1, (85) 

r=l 

see Eq. ( 68 ). Accordingly, the partition Eq. (70) becomes 

[Nt] = AtBtDt... A2B2D2A1B1D1AQ. ( 86 ) 

We shall use terms A-qubit, 5-qubit, and 5-qubit to indicate which of the blocks in Eq. ( 86 ) contains 
a given qubit. Note that the site w^®’ of Dr and the site of Dr-i coincide, see Fig. 6 . However, 
since we placed only one qubit at and w^, the total number of qubits per site is at most two. 

For each r = 1,..., t dehne 2r additional gauge generators gl and h* with i = 1,..., r as shown 
below, see also Fig. 8 . 

gl = {w^^ + w^^^^)[Dr], (87) 

where i = 1,... ,r — 1, 

gr = {w^ + w^'')[Dr], ( 88 ) 
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(89) 


K = (m”-‘ + «i“)|£>J + + u^')[Br] + (k“-‘ + K“)l'4r-i], 

where i = 1,..., r — 1, and 

h; = ^ (^2r-l ^ ^2r ^ u^r+l~^[Br\ + [/l,_i]. (90) 


9! 



Figure 8: Additional gauge generators h\. and gl for r = 3. The generator h* acts on i?-qubits at 
the boundary of (blue) and A-qubits at the boundary of Ar_i (red). In addition, h® acts on two 
D-qubits (hlled circles). The generator gl acts only on two D-qubits. The special faces 6® and c® of 
the color code lattice are used to dehne modihed stabilizers. The original stabilizer U[Ar\ + U[Br\ G % 
must be replaced by gl + 6 ®[Ar] + &®[i?r] G Ut- The original stabilizer B[Ar] + B[Br] G Tt must be 
replaced by g'!;._^^ + B[Ar] + c®[i?r] G Ut- 

Note that all additional generators except for g^ are spatially local. Furthermore, 

r 

^r,r-l = h® + g).. (91) 

i=l 

This identity will allow us to get rid of the non-local generators Ur^r-i- More formally, dehne a 
subspace Ut C such that 


t 

Ut = % + Y,{9l---^9iyr.---.K)- 

r=l 


(92) 


Note that Ut uniquely dehnes Ut since applying the dot operation twice gives the original subspace, 
see Section 3. We shall see that Ut can be regarded as an extended version of 7*. Similarly, dehne a 
subspace Vt C such that 


t 

Vt = Ct + '^{9l, ■ ■ ■ ,9l,hl,..., K). 

r=l 


(93) 
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We shall see that Vt can be regarded as an extended version of Ct. Here it is nnderstood that vectors 
from Tt and Ct are extended to D-qnbits by zeros. From Eq. (82) we infer that 

Vt C Ut and Ut^Vt. (94) 

Let ns establish some basic properties of Vt and Ut. 

Lemma 5. Ut ^ Ut- Furthermore, if h ^ Ut then the restriction of h onto the union of A and B 
qubits is contained in 7*. 

Lemma 6. Vt F Vt. Furthermore, if h ^ Vt then the restriction of h onto the union of A and B 
qubits is contained in Ct. 

Since the proof of the two lemmas is identical, we only prove Lemma 5. 

Proof. Snppose h EUt. By dehnition of the dot operation, h has even weight and orthogonal to any 
vector of Ut. Let hn be a restriction of h onto the nnion of all D-qnbits. Since hjj is orthogonal 
to any generator g^, we infer that ho is a. linear combination of gl. This implies h^ E Ut. Let 
be a restriction of h onto the nnion of all A- and H-qnbits. Since h E S and hr, E 8, we infer that 
hAB ^ £. The inclusion Tt F Ut implies Hab £ q H T = 7*. This proves the second statement. 
Finally, Tt EL Tt ^ Ut implies hAB E Ut, that is, h E Ut. Tl 

Combining the lemmas and Eq. (94) yields 

UtFVtCVtC Ut. (95) 

One can view Eq. (95) as an extended version of Eq. (82). We can now dehne extended versions of 
the doubled color codes, see Table 3. 



Transversal gates 

8tabilizer group 

Gauge group 

C-code 

Clihord group 

088 ( 7 ) 4 , 7 ) 4 ) 

088 ( 7 ) 4 , 7 ) 4 ) 

T-code 

T gate 

088(^4, W 4 ) 

088(^4,4) 

Base code 


088 (^ 4 , 7 ) 4 ) 

088 ( 7 ) 4 ,^ 4 ) 


Table 3: Extended doubled color codes. 


As before, we dehne the base code such that its stabilizer group is the intersection of stabilizer 
groups of all codes in the family. Lemma 6 implies that Vt is self-orthogonal, that is, the group 
088 ( 7 ) 4 , 7 ) 4 ) is abelian and the C-code is well-dehned. Likewise, Eq. (95) implies that Ut and Vt 
are mutually orthogonal, so that the group CSS(Ut,Vt) is abelian and the base code is well-dehned. 
8ince we already know that Ct is doubly even and T is triply even. Lemmas 5,6 have the following 
corollary. 
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Corollary 3. The subspace Ut is triply even with respect to the same subsets as %. The subspace Vt 
is doubly even with respect to the same subsets as Ct- 

This shows that the extended C and T codes have transversal Clifford gates and the T-gate 
respectively. Let us describe gauge generators of the extended C-code. From Eq. (91) we infer that 
the non-local generators can be removed from the generating set of Vt. Thus the extended 

C-code has only face-type gauge generators and the additional generators gl,..., g'!^, hi,..., h'^. The 
latter are spatially local except for g^.. A direct inspection of Eq. (93) shows that generators of Vt 
can be partitioned into mutually disjoint subsets supported on regions 

Alt = At, Mt-i = BtDtAt-i, Mt-2 = Bt-iDt-iAt-2, ... , Mq = B\DiAq. 

Thus we can analyze properties of the extended C-code on each region separately. Generators of 
Dt supported on Mt have a form f[At\ with f & St. These generators describe the regular color code 
St- Thus the extended C-code can be converted to the regular color code CSS{St,St) by discarding 
all the regions except for M*. As before, the extended T-code is obtained from the extended C-code 
by adding edge-type Z-stabilizers. 

Next let us prove that the extended codes have distance 2t + 1. Combining Eq. (7) and Eq. (95) 
one can see that the distance of any code defined in Table 3 is lower bounded by d{Ut)- Thus it 
suffices to prove that diUt) = 2t + 1. We shall need an explicit expression for generators of Ut. Note 
that 

r 

9r = ^ Alt. (96) 

^=1 

Indeed, gr does not overlap with generators of Tt and has even overlap with all additional generators 
gl and h*. Next consider a generator gl with 1 < i < r. Let 6* G Sr be the face of located directly 
below gl, see Fig. 8. We claim that 

131 = gl + -|- V[Br] G Ut (97) 

for all r = 1,... ,t and alH = 1,..., r — 1. Indeed, /3* has even overlap with all generators of % since 
6*[Ar] -|- b^[Br] G 7*. Furthermore, (31 has even overlap with all additional generators g].,. It remains 
to check that (31 has even overlap with the additional generators h*,. The only non-trivial case is h^ 
with j = ioij = i + l, see Fig. 8. In this case both gl and h^[Ar] -|- B[Br\ have odd overlap with hi, 
so that (31 has even overlap with h^. This proves that G W*. 

Likewise, let c* G St-i be the face of Aj,_i located directly above gl, see Fig. 8. The same 
arguments as above show that 


ll = qIA d[Ar-i] d[Br-\\ G Ut (98) 

for all r = 2,..., t and for alH = 1,..., r — 1. Note that the sublattice At has only special faces 6* 
whereas Ai has only special faces d. All other sublattices A^ have both types of special faces. Now 
we are ready to describe generators of Ut. 
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Lemma 7. Suppose 2 < r < t — 1. Let S* C Sr be the subspace spanned by all faces of the color 
code lattice Ar except for the special faces 6* and c\ Let S^ C St be the subspace spanned by all faces 
of At except for b\ Let SI C Si be the subspace spanned by all faces of Ai except for c^. Then 

t 

Ut = + f[Br] : / e 5;) + {BrAr-l) 

r=l 

t 

+ ^(s'r) + {fdl, 7r • * = 1, • • •, ?" — 1). (99) 

r=l 

Proof. We have already shown that the last two terms in Eq. (99) are contained in Ut- A direct 
inspection shows that vectors f[Ar] + f[Br] with f & S* and BrAr-i have even overlap with all 
additional generators g^, h*. Furthermore, f[Ar] + f[Br] and BrAr-i are contained in Tt and thus 
have even overlap with any vector in 7*. This proves the inclusion 3 in Eq. (99). 

Conversely, consider any k ^Ut- We have to prove that k is contained in the righthand side of 
Eq. (99). The same arguments as in the proof of Lemma 5 show that 

t r 

k = kAB + kn, kAB ^ ko = EE xlgi, (100) 

r=l i=l 

where and ko have support only on Ai^-qubits and H-qubits respectively. Here x* G F 2 are 
some coefficients. Let us modify k according to 

t r 

k ^ k + j2E^M- 

r=l i=l 

We still have the inclusion k EUk since ElLk, see above. Furthermore, the term x*/?* cancels the 
term x^gl in kn and modihes the term kAB according to kAB t— kAB + h^[Ar\ + P[Br\. Now k has 
support only on HH-qubits and k E Tt- Let us express fc as a sum of generators of Tt defined in 
Eq. (71). This yields 

t r—1 

k = k* + y ^ yl{B[Ar] + b^[Br]) + zl{c^[Ar-l] + C*[i?,._i]) 

r=2 i=l 

for some 

t 

k* E '^^{f[Ar] + f[Br] : f E S*) + (BrAr-l) 

r=l 

and some coefficients yl,zl E F 2 . Note that k* is contained in the righthand side of Eq. (99). 
Furthermore, k* G Ut which implies k + k* E Ut- In particular, k + k* has even overlap with all 
generators hi- On the other hand, hi has odd overlap with P[Ar\ + P[Br\ and with P[Ar-i\ +P[Br-i] 
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for i = j, j — 1, see Fig. 8 . Thus k + k* may have even overlap with hi only if yl = z* for all i and r. 
Using the identity 

+ i>^[Br\) + (c*[^r-l] + C*[i?r--l]) = /^r + 7r 

we conclude that k + k* is contained in the last term in Eq. (99). Since k* is contained in the sum 
of the hrst two terms in Eq. (99), we have proved the inclusion C in Eq. (99). □ 

The following lemma is the most difficult part of the proof. 

Lemma 8. d{Ut) = d{Tt) = 2t + 1. 

Proof. We shall use induction in t. The base of induction is t = 1. In this case there are only two 
D-qubits, Di = and one additional gauge generator gi = g\ = + w^. Furthermore, gi 

is a stabilizer, gi G Wi. Consider a minimum weight vector k G Uf fl O such that d{Ui) = \k\. Since 
k has even overlap with ^fi, one has either gi k oi giP k = The hrst case can be ruled out 
since k + gi & Uf fl O would have weight \k\ — 2. Thus k has support only on Ai^-qubits. Lemma 5 
implies that Ui and 7i have the same restriction on Ai^-qubits, that is, d{Ui) = d(7i) = 3. 

Consider now some t > 2. Let us rewrite the partition in Eq. (86) as 

[W] = AfBtDtCt, Ct = At-iBt-iDt-i ... AiBiDiAq. 

Consider an arbitrary vector k G Uf- fl O. Let us write 

k = a[At\ + l3[Bi\ + d[Di\ + 

for some vectors 

a,/3GF™‘, ^gF^*, 7 GFf-L 

First we claim that 

a G (9, /3 + 7 G F, and 5 G F. 

Indeed, since k has even overlap with the stabilizer gt with supp( 5 fi) = Dt, see Eq 
0 = k'^gt = |5| (mod 2), that is, 6 & S. Furthermore, Lemma 7 implies 

t t-i 

h = BtCt = BrAr-l + gr G Ut- 

r=l r=l 

Since supp(h) = BtCt and 0 = kdh = |/?| + | 7 | (mod 2), we infer that /? + 7 G F. Finally, a G (9 
follows from the above and the assumption that k E O. 

We shall refer to a substitution k ^ k + g with g E Ut as a gauge transformation. Note that 
gauge transformations preserve the set Uf- fl O. Our strategy will be to choose a sequence of gauge 
transformations that transform k into a form 

k = P[Af\ + P[Bf\ + (102) 


( 101 ) 

. (96), we infer that 
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without increasing the weight of k. Here e® is some basis vector of F™*. Let us hrst assume that k 
has the desired form Eq. (102) and show that this implies \k\ >2t + l. Indeed, using Eq. (99) one 
can check that e®[Hj] + e®[i?t] is orthogonal to all generators of Ut except for BtAt-i. Therefore if 
k G has a special form Eq. (102) then 7 G U^_i- Furthermore, k ^ O implies 7 G (9, that is, 

7 G n O. By induction hypothesis, I 7 I > 2f — 1 and thus \k\ > 2t + 1. 

It remains to show that k can be transformed into the desired form Eq. (102) without increasing 
the weight. First, choose any i G supp(Q!) and let a' = a + e\ Note that a' E £ due to Eq. (101), so 
that a'[Ai\ + «'[!?*] ElAt- In addition, |q!'| = |a| — 1, so that 

\k + = 1 + |/? + + I7I + |( 5 | < 1 + \a \ + |/?| + I7I + 1*^1 = + |/ 3 | + Ibl + I'^l = 1^1- 

Thus we can transform k io k + + a'[Bi\ obtaining 

k = + /3[i?t] + S[Dt\ + 'y[Ct] 

for some new vector f3 G F™*. Dehne 6 * = /3 + e® so that 

k = e®[Hj] + e®[Ht] + 0[B^\ + 5[D^\ + (103) 

Consider hrst the case when 0 = 0. Since k has even overlap with the stabilizer jSl and so does 
e®[Ht] + e®[i?t], we conclude that S[Di\ has even overlap with the gauge generator gl for any j. This 
is possible only if 5 = 0 or 5 = 1. If 5 = 0 then k already has the desired form Eq. (102). If 5 = Iw 
then k + gl has the desired form and \k + gl\ < \k\. 

Consider now the case 0 7 ^ 0. Since k has even overlap with the stabilizers /Sj and so does 
e®[y4t] +e®[Ht], we conclude that S[Dt\ has even (odd) overlap with a generator gl iff 0 has even (odd) 
overlap with the special face V, see Fig. 8 . Let e(0) = |0| (mod 2) be the total parity of 0. Let a{6) 
be the set of faces of that have odd overlap with 0. The above shows that 

ct(0) C {61,62^..., (104) 

Note that 

|a(0)|=O (mod 2) (105) 

since 6 E £, see Eq. (101). We claim that one can choose 6* C At such that 

1^*1 < l^'l, a{6*) = a{6), e( 0 *) = e( 0 ), and supp(0*) C {n^, ^ (106) 

Recall that u^,u^,... are the sites of A* located on the boundary facing Ai_i, see Fig. 7. 

Indeed, it is straightforward to choose some 6* that satishes all above conditions except for the hrst 
one. However, the property that 9* is supported on the boundary oj = implies 

that its weight cannot be decreased by adding faces of the lattice At. Indeed, the boundary oj is 
a minimum weight logical operator of the color code. This implies |a; + /| > |a;| for any stabilizer 
f E St■ Equivalently, |/\ci;| > |/na;|. This implies 

\f + e*\ = \f\uj\ + \{fnu) + e*\ > \fnu\ + \ifnu) + e*\ > \e*\. 
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Finally, we note that cr{6*) and e{6*) fix 6* modulo stabilizers / G St- This proves that we can satisfy 
all conditions in Eq. (106). 

Note that 0[Bt\ + 0*[Bt\ G Ut since 9 + 9* G St- A gauge transformation k ^ k + 9[Bt\ + 9*[Bt\ 
can potentially increase the weight of A; if i G 0 but i ^ 9*. However, the weight can increase at most 
by two. We will “borrow” two units of weight keeping in mind that at least one of the subsequent 
gauge transformations has to decrease the weight of k at least by two, so that we maintain a zero 
weight balance. Transforming k to k + 9[Bt] + 9*[Bt] we obtain 

k = e*[At] + e*[i?t] + 9*[Bi\ + 5[Dt\ + (107) 

We have to two consider two cases. 

Case 1: 9* G £. Then ^ 9* and ^ 9*, so that 9* consists of disjoint paths connecting 

consecutive pairs of faces in (t{9*). Consider any such path 

TT = + „2i+2 + _ _ _ + „2m-l ^ ^2m ^ q* 

that creates a pair of syndromes at faces 6® and 6"® for some i < m. An example of such path with 
i = 1 and m = 4 is shown on Fig. 9. As we argued above, d[Di\ must have odd overlap with gl and 
g^. Applying, if necessary, gauge transformations k ^ k + gl and fc —)• /c + we can assume that 
y^ 2 i+i g ^ ^ 2 m g Then a gauge transformation 

m m —1 

k ^ k + ^ ^ gl 

p=i+l p=i-\-l 

cleans k out of all qubits vr[i?i], qubits + tc^”®)[Zli], and potentially adds weight at |7r| qubits 

(^2i+l ^2i+2 ^ .j;2m-l _j_ 

Overall, the weight decreases at least by two, see Fig. 9. 

Case 2: 9* G O. Then G 9* and G 6**, so that 9* consists of a path connecting the leftmost 
face in a(6**) to the site a path connecting the rightmost face in o'{9*) to the site u^, and, 

possibly, disjoint paths connecting consecutive pairs of faces in (j{9*), see Fig. 10. The latter can 
be cleaned out in the same fashion as in Case (1), so below we focus on the former. Let 6* be the 
rightmost face in o'{9*). Then 9* contains a path 

'bright = + . . . + ^ 

Let 6™ be the leftmost face in 'y{9*). Then 9* contains a path 

= ^2^+1 ^ ^ ^ ^^2t ^ ^24+1 _ 

As we argued above, S[Dt] must have odd overlap with gl and g^. Applying, if necessary, gauge 
transformations k ^ k + gl and k ^ k + g"^ we can assume that G <5 and i(; 2 m+i ^ Then a 
gauge transformation 

k ^ k + ^ hi + ^ hi + ^ gl + ^ gl 

p=l p=m-\-l p=l p=m-\-l 


47 



Figure 9: Example of a vector k defined in Eq. (107) such that 6 * E S creates a pair of syndromes 
at faces and b^. Sites of 6 * are indicated by red circles. In order to have zero syndrome for 
stabilizers and the support of k must also include sites and (possibly, after a gauge 
transformation k ^ k + g] and k ^ k + gf)- The sites and are indicated by blue circles. A 
gauge transformation k ^ k + hf + + hf + g^ + g^ cleans out k from all qubits indicated by blue 

and red circles, potentially adding weight to qubits indicated by green circles. Overall, the weight of 
k decreases at least by two. Here t = 5. 

cleans k out of all qubits of {irieft + '^right)[Bt\, qubits (tc^* and potentially adds weight 

at \7Tright\ qubits (n^ + ... + and \7rieft\ — 2 qubits (n^m+i Overall, the 

weight decreases at least by four, see Fig. 10. In both cases, we transform k to the desired form and 
the weight decreases at least by two. □ 

The hnal weight reduction step is to transform the long-range generators g^. into a local form. 
We shall use a coding theory analogue of the subdivision gadget used to simulate long-range spin 
interactions by short-range ones [26] . Consider an arbitrary subspace U such that U <El{ and n 
is odd. It defines a CSS code with a gauge group CSS(W,W) and a stabilizer group CSS(7/,7/). Note 
that 7/ C 7/ implies that 77 is self-orthogonal, so that the stabilizer group is abelian and the code is 
well-defined. Suppose 77 contains a weight-two vector supported on the first two qubits, -t- G 77. 
We envision a scenario when qubits 1 and 2 occupy two remote lattice locations (for example sites 
and in the above construction), such that -f- is not spatially local. Let us add two ancillary 
qubits labeled a and b. We envision that all vectors -|- e“, e“ -1- e^, and are spatially local 

(or, at least, more local compared with -f e^). Define a subspace V C 82 '^“^ such that 

V= (el + e“,e“ + e^e^ + e^)+77. (108) 

Here it is understood that vectors of 77 are extended to the ancillary qubits by zeroes. Using analogues 
of Lemmas 5,6 one can easily show that V is self-orthogonal, so that V defines an abelian stabilizer 
group CSS(V, V). Below we prove that the codes CSS(77,77) and CSS(V, V) have the same distance 
and the same transversality properties. 
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Figure 10: Example of a vector k defined in Eq. (107) such that 6 * & O creates a pair of syndromes 
at faces and b^. Sites of 6* are indicated by red circles. In order to have zero syndrome for 
stabilizers (3^ and the support of k must also include sites and (possibly, after a gauge 

transformation k ^ k + gf and k ^ k + gl)- The sites and w'^ are indicated by blue circles. A 
gauge transformation k ^ k + h\ + + g] + gj + gl cleans out k from all qubits indicated 

by blue and red circles, potentially adding weight to qubits indicated by green circles. Overall, the 
weight of k decreases at least by four. Here t = 5. 

Lemma 9 (Subdivision gadget). Consider the subspaces U and V as above. Suppose d{U) > 3. 
Then d{V) = d{U). Furthermore, ifU is triply (doubly) even with respect to some subsets C [n] 
then V is triply (doubly) even with respect to the same subsets. 

Proof. Let U 12 and U 12 be the subspaces including all vectors g eU and g eU respectively such that 
91 = 92 = 0, that is, g does not include qubits 1 and 2. By assumption, G W and thus /i = /2 

for any / G 7/ since / must have even overlap with any element of U. There must be at least one 
vector f eU such that /i = /2 = 1 since otherwise d(U) = 1. Thus U can be represented as 

U = {e^ + + g) + U 12 for some g G 7/i2. (109) 

We claim that 

V = (h) + U 12 , where h = + e“ + + 5 f. (HO) 

Let us hrst prove the inclusion D in Eq. (HO). We have h G since + g E U C and 

e“ + G since no vector in lA includes a or b. Taking into account Eq. (108) one gets h E V-*- 
and thus h E V. The inclusion V(i 2 C V follows trivially from the dehnitions. Next let us prove the 
inclusion C in Eq. (HO). Consider any vector k E V. We note that ki = ka = kb = /c 2 since k must 
have even overlap with any vector in V. Replacing, if necessary, k hy k + h we can assume that 
ki = ka = kb = k 2 = 0. Since W C V we infer that k E lA^. Taking into account that k has even 
weight one gets k ElA = LA, that is, k G IA 12 . We have proved Eq. (HO). 
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Suppose / G V-*- n C> is a minimum weight vector such that |/| = diV). Then / includes at most 
one of the qubits 1, a, 6, 2 since otherwise we would be able to reduce the weight of / by a gauge 
transformation / ^ / + /^ where /' is contained in the hrst term in Eq. (108). Without loss of 
generality fa = fb = f 2 = 0. From Eq. (110) we infer that h'^f = 0 and / G Since fa = ft = 0, 
this implies that / is orthogonal to + 5 ^ and thus / G see Eq. (109). This shows that 

/ G lA^ n O, that is, d{V) = |/| > d{U). The opposite inequality, diV) < d{U), is obvious since 
extending any vector / G fl by zeroes to qubits a and b gives a vector f E H O. The last 
statement of the lemma follows from the fact that U and V have the same restriction on any subset 
C [u], see Eq. (110). □ 

It remains to apply the subdivision gadget to the long-range gauge generators = {w^ +w‘^'^)[Dr\, 
where r = 2,... ,t (note that g\ is already spatially local). Let us add a second qubit at each site tc* 
of the region D^, except for and We shall denote these extra qubits as w^,... such 

that Dr = {w^yw"^,... Qubits w\w^ share the same site of the lattice. The 

total number of qubits becomes 

t 

Kt = Nt + 5^(2r - 2) = 2t^ + + 6t + 1. 

r=2 


Dehne a subspace Vt C such that 

t 2r-2 

= Dt + + w‘^)[Dr]) + ^ + w'^^)[Dr]) + + W^~^^)[Dr]). (HI) 

r=2 i=2 

Here it is understood that vectors of Ut are extended to the added qubits by zeroes. Note that all 
generators of Vt are spatially local since the long-range generator gl can be decomposed as 

gl = + w'^'' = {wi -f W 2 ) + {w2 + W 3 ) + • • • + {W2r-1 + W2r), (112) 

where each term has support on a single face of the lattice. Here we omitted [Dr] to simplify notations. 
We shall regard Vt as an extended version of Ut. Similarly, dehne a subspace Dt Q £ 2 ^ such that 

t 2r-2 

Dt = 'Dt + +iv‘^)[Dr]) + + w‘^'^)[Dr]) + ^((w* (113) 

r=2 i=2 

We shall regard Dt as an extended version of Dt. Note that all generators of Dt are spatially local. 
Using Eq. (95) and analogues of Lemmas 5,6 one can easily show that 

W C C C Vt. (114) 

We are now ready to dehne the hnal version of the doubled color codes, see Table 4, that satisfy all 
the properties announced in Section 8 (up to relabeling of the subspaces Ct, % into Vt). Combining 
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Eq. (7) and Eq. (114) one can see that the distance of any code dehned in Table 4 is lower bounded 
by diVt). A recursive application of the subdivision gadget lemma shows that diVt) = d{Ut) = 2t +1. 
Furthermore, Vt is triply-even with respect to the same subsets as Ut and is doubly even with 
respect to the same subsets as T>t. Thus the hnal C and T codes have transversal Clifford gates 
and the T-gate respectively. Let us describe gauge generators of the hnal C-code. From Eq. (112) 



Transversal gates 

Stabilizer group 

Gauge group 

C-code 

Clihord group 


css{A,A) 

T-code 

T gate 

CSS(W,W) 

CSS(W,W) 

Base code 


CSS{Vt,J^t) 

CSS(T-i,W) 


Table 4: Final version of the doubled color codes. All gauge generators of the C-code are spatially 
local. The T-code can be obtained from the C-code by measuring syndromes of edge-type stabilizers. 
The base code is dehned such that its stabilizer group is the intersection of all other stabilizer groups. 

we infer that the long-range generators g^. can be removed from the generating set. Thus the hnal 
C-code has only face-type gauge generators, the additional generators gl,..., gr~^, hi,, h^, and 
the additional generators that appear in Eq. (113). All these generators are spatially local. The same 
arguments as above show that a restriction of the hnal C-code onto the region At coincides with the 
regular color code CSS{St,St). As before, the hnal T-code is obtained from the extended C-code by 
adding edge-type stabilizer generators. 

Finally, let us point out that the gauge generator € 7t is already spatially local, so that we 
do not have to apply the weight reduction steps to cui^o- Thus the total number of physical qubits 
can be reduced from Kt to Kt — 2 = 2t^ + 8t‘^ + 6t — 1. 


Acknowledgments 

SB thanks Andrew Landahl for fruitful discussions and helpful suggestions at the early stages of this 
project. SB acknowledges NSF grant PHY-1415461. 


References 

[1] R. Barends, J. Kelly, A. Megrant, A. Veitia, D. Sank, E. Jehrey, T.C. White, J. Mutus, A. Fowler, 
B. Campbell, et ah Superconducting quantum circuits at the surface code threshold for fault 
tolerance. Nature, 508(7497):500-503, 2014. 

[2] J. Kelly, R. Barends, A. Fowler, A. Megrant, E. Jehrey, T.C. White, D. Sank, J. Mutus, 
B. Campbell, Y. Chen, et al. State preservation by repetitive error detection in a supercon¬ 
ducting quantum circuit. Nature, 519(7541):66-69, 2015. 


51 










[3] A. Corcoles, E. Magesan, S. Srinivasan, A. Cross, M. Steffen, J. Gambetta, and J. Chow. 
Demonstration of a quantum error detection code using a square lattice of four superconducting 
qubits. Nature Communications, 6:6979, 2015. 

[4] E. Dennis, A. Kitaev, A. Landahl, and J. Preskill. Topological quantum memory. J. of Math. 
Phys., 43(9):4452-4505, 2002. 

[5] R. Raussendorf and J. Harrington. Fault-tolerant quantum computation with high threshold in 
two dimensions. Phys. Rev. Lett., 98(19):190504, 2007. 

[6] A. Fowler, A. Stephens, and P. Groszkowski. High-threshold universal quantum computation 
on the surface code. Phys. Rev. A, 80(5):052312, 2009. 

[7] C. Jones. Logic synthesis for fault-tolerant quantum computers, preprint arXiv:1310.7290, 2013. 

[8] A. Fowler, S. Devitt, and C. Jones. Surface code implementation of block code state distillation. 
Scientific Reports, 3:1939, 2013. 

[9] B. Eastin and E. Knill. Restrictions on transversal encoded quantum gate sets. Phys. Rev. Lett., 
102(11):110502, 2009. 

[10] S. Bravyi and R. Konig. Classihcation of topologically protected gates for local stabilizer codes. 
Phys. Rev. Lett, 110(17):170503, 2013. 

[11] F. Pastawski and B. Yoshida. Fault-tolerant logical gates in quantum error-correcting codes. 
Phys. Rev. A, 91(1):012305, 2015. 

[12] M. Beverland, R. Konig, F. Pastawski, J. Preskill, and S. Sijher. Protected gates for topological 
quantum held theories, preprint arXiv:1409.3898, 2014. 

[13] A. Paetznick and B. Reichardt. Universal fault-tolerant quantum computation with only 
transversal gates and error correction. Phys. Rev. Lett., 111(9):090505, 2013. 

[14] J. Anderson, G. Duclos-Cianci, and D. Poulin. Fault-tolerant conversion between the Steane 
and Reed-Muller quantum codes. Phys. Rev. Lett., 113(8):080501, 2014. 

[15] E. Knill, D. Leibfried, R. Reichle, J. Britton, R.B. Blakestad, J.D. Jost, C. Langer, R. Ozeri, 
S. Seidelin, and D.J. Wineland. Randomized benchmarking of quantum gates. Phys. Rev. A, 
77(1):012307, 2008. 

[16] J. Chow, J. Gambetta, L. Tornberg, J. Koch, L. Bishop, A. Houck, B.R. Johnson, L. Frunzio, 
S. Girvin, and R. Schoelkopf. Randomized benchmarking and process tomography for gate 
errors in a solid-state qubit. Phys. Rev. Lett., 102(9):090502, 2009. 

[17] B.J. Brown, H. N. Naomi, and D.E. Browne. Fault-tolerant error correction with the gauge 
color code, preprint arXiv:1503.08217, 2015. 


52 



[18] H. Bombin. Gauge color codes: optimal transversal gates and gauge fixing in topological stabi¬ 
lizer codes. New J. Phys., 17(8):083002, 2015. 

[19] C. Wang, J. Harrington, and J. Preskill. Confinement-Higgs transition in a disordered gauge 
theory and the accuracy threshold for quantum memory. Ann. of Phys., 303(l):31-58, 2003. 

[20] K. Betsumiya and A. Munemasa. On triply even binary codes, preprint arXiv:1012.4i34, 2010. 

[21] A. R. Calderbank and P. W. Shor. Good quantum error-correcting codes exist. Phys. Rev. A, 
54(2):1098, 1996. 

[22] A.M. Steane. Multiple particle interference and quantum error correction. Proc. Roy. Soc. Lond. 
A, 452:2551-2577, 1996. 

[23] H. Bombin and M. A. Mart in-Delgado. Topological quantum distillation. Phys. Rev. Lett., 
97(18):180501, 2006. 

[24] A. Kubica and M. Beverland. Universal transversal gates with color codes: a simplified approach. 
Phys. Rev. A, 91(3):032330, 2015. 

[25] A. Kitaev. Fault-tolerant quantum computation by anyons. Ann. of Phys., 303(l):2-30, 2003. 

[26] R. Oliveira and B. Terhal. The complexity of quantum spin systems on a two-dimensional square 
lattice. Quant. Inf. and Comp., 8(10):900-924, 2008. 

[27] S. Bravyi and J. Haah. Magic-state distillation with low overhead. Phys. Rev. A, 86(5):052329, 
2012 . 

[28] A. Landahl and G. Ryan-Anderson. Quantum computing by color-code lattice surgery, preprint 
arXiv:1407.5103, 2014. 

[29] V. Kliuchnikov, D. Maslov, and M. Mosca. Fast and efficient exact synthesis of single-qubit 
unitaries generated by Glifford and T gates. Quant. Inf. and Comp., 13(7-8):607-630, 2013. 

[30] P. Selinger. Efficient Glifford-|-T approximation of single-qubit operators. Quant. Inf. and 
Comp., 15(1-2):159-180, 2015. 


53 



